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SIGNAL PROCESSING APPARATUS AND METHOD 

Field of the Invention 
The present invention relates to the field of signal 
processing. More specifically, the present invention relates 
5 to the processing of measured signals to remove undesired 

portions when little is known about the undesired signal 
portion. 

Background of the Invention 
Signal processors are typically employed to remove 

10 undesired portions from a composite measured signal including 
a desired signal portion and an undesired signal portion. If 
the undesired signal portion occupies a different frequency 
spectrum than the desired signal* .then conventional filtering 
techniques such as low pass, band pass, and high pass 

15 filtering could be used to separate the desired portion from 

the total signal. Fixed single or multiple notch filters 
could also be employed if the undesired signal portions (s) 
exist at a fixed frequency (s) . 

However, it is often the case that an overlap in 

20 frequency spectrum between the desired and undesired signal 

portions does exist and the statistical properties of both 
signal portions change with time. In such cases, 
conventional filtering techniques are totally ineffective in 
extracting the desired signal. If, however, a description of 

25 the undesired portion can be made available,' adaptive noise 
canceling can be employed to remove the undesired portion of 
the signal leaving the desired portion available for 
measurement. Adaptive noise cancelers dynamically change 
their transfer function to adapt to and remove the undesired 

30 signal portions of a composite signal. Adaptive noise 
cancelers require a noise reference signal which is 
correlated to the undesired signal portion. The noise 
reference signal is not necessarily a representation of the 
undesired signal portion, but has a frequency spec trum which 
35 is similar to that of the undesired signal, in many cases, 




it requirW^ considerable ingenuity ^fe determine a noise 
reference signal since nothing is a priori known about the 
undesired signal portion. 

One area where composite measured signals comprise a 
5 desired signal portion and an undesired signal portion about 

which no information can easily be determined is 
physiological monitoring. Physiological - monitoring 

apparatuses generally measure signals derived from a 
physiological system, such as the human body. Measurements 
10 which are typically taken with physiological monitoring 
systems include electron cardiographs, blood pressure, blood 
gas saturation (such as oxygen saturation), capnographs, 
heart rate, respiration rate, and depth of anesthesia, for 
example. Other types of measurements include those which 
15 measure the pressure and quantity of a substance within the 
body such as breathalizer testing, drug testing, cholesterol 
testing, glucose testing, arterial carbon dioxide testing, 
protein testing, and carbon monoxide testing, for example. 
The source of the undesired signal portion in these 
20 measurements is often due to motion of the patient, both 
external and internal (muscle movement, for example) , during 
the measurement process. 

Knowledge of physiological systems, such as the amount 
of oxygen in a patient's blood, can be critical, for example 
25 during surgery. Data can be determined by a lengthy invasive 
procedure of extracting and testing matter, such as blood, 
from a patient, or by more expedient, non-invasive measures. 
Many types of non-invasive measurements can be made by using 
the known properties of energy attenuation as a selected form 
30 of energy passes through a medium. 

Energy is caused to be incident on a medium either 
derived from or contained within a patient and the amplitude 
of transmitted or reflected energy is then measured. The 
amount of attenuation of the incident energy caused by the 
35 medium is strongly dependent on the thickness and composition 
of the medium through which the energy must pass as well as 
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the specific form of energy selected Information about a 
physiological system can be derived from data taken from the' 
attenuated signal of the incident energy transmitted through 
the medium if the noise can be removed. However, non- 
invasive measurements often do not afford the opportunity to 
selectively observe the interference causing the undesired 
signal portion, making it difficult to remove. 

These undesired signal portions often originate from 
both AC and DC sources. The first undesired portion is an 
easily removed DC component caused by transmission of the 
energy through differing media which are of relatively 
constant thickness within the body, such as bone, tissue, 
skin, blood, etc. Second, is an erratic AC component caused 
when differing media being measured are perturbed and thus, 
change in thickness while the measurement is being made. 
Since most materials in and derived from the body are easily 
compressed, the thickness of such matter changes if the 
patient moves during a non-invasive physiological 
measurement. Patient movement can cause the properties of 
energy attenuation to vary erratically. Traditional signal 
filtering techniques are frequently totally ineffective and 
grossly deficient in removing these motion induced effects 
from a signal. The erratic or unpredictable nature of motion 
induced undesired signal components is the major obstacle in 
removing them. Thus, presently available physiological 
monitors generally become totally inoperative during time 
periods when the patient moves. 

A blood gas monitor is one example of a physiological 
monitoring system which is based upon the measurement of 
energy attenuated by biological tissues or substances . Blood 
gas monitors transmit light into the tissue and measure the 
attenuation of the light as a function of time. The output 
signal of a blood gas monitor which is sensitive to the 
arterial blood flow contains a component which is a waveform 
representative of the patient's arterial pulse. This type of 
signal, which contains a component related to the patient's 
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.^^alled a plethysmographi^^L 



pulse, i^^alled a plethysmographic^/ave, and is shown in 
Figure 1 as curve Y. Plethysmography waveforms are used in" 
blood pressure or blood gas saturation measurements, for 
example. As the heart beats the amount of blood in the 
5 arteries increases and decreases, causing increases and 

decreases in energy attenuation, illustrated by the cyclic 
wave Y in Figure- 1. 

Typically, a digit such as a finger, an ear lobe, or 
other portion of the body where blood flows close to the 
10 skin, is employed as the medium through which light energy is 
transmitted for blood gas attenuation measurements* The 
finger comprises skin, fat, bone, muscle , etc., shown 
schematically in Figure 2, each of which attenuates energy 
incident on the finger in a generally predictable and 
15 constant manner. However, when fleshy portions of the finger 
are compressed erratically, for example by motion of the 
finger, energy attenuation becomes erratic. 

An example of a more realistic measured waveform S is 
shown in Figure 3, illustrating the effect of motion. The 
20 desired portion of the signal Y is the waveform 
representative of the pulse, corresponding to the sawtooth- 
like pattern wave in Figure 1. The large, motion-induced 
excursions in signal amplitude hide the desired signal Y. : It 
is easy to see how even small variations in amplitude make it 
25 difficult to distinguish the desired signal Y in the presence 
of a noise component n. 

A specific example of a blood gas monitoring apparatus 
is a pulse oximeter which measures the saturation of oxygen 
in the blood. The pumping of the heart forces freshly 
30 oxygenated blood into the arteries causing greater energy 
attenuation. The saturation of oxygenated blood may be 
determined from the depth of the valleys relative to the 
peaks of two plethysmography waveforms measured at separate 
wavelengths. However, motion induced undesired signal 
35 portions, or motion artifacts, must be removed from the 
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measured ^^nal for the oximeter to continue the measurement 
during periods when the patient moves. 

Summary of the Invention 
The present invention is a signal processor which 
5 acquires a first signal and a second signal that is 

correlated to the first signal. The first signal comprises 
-a first desired signal portion and a first undesired signal 
portion. The second signal comprises a second desired signal 
portion and a second undesired signal portion. The signals 

10 may be acquired by propagating energy through a medium and 
measuring an attenuated signal after transmission or 
reflection. Alternatively, the signal may be acquired by 
measuring energy generated by the medium. 

The first and second measured signals are processed to 

15 generate a noise reference signal which does not contain the 
desired signal portions from either of the first or second 
measured signals. The remaining undesired signal portions 
from the first and second measured signals are combined to 
form a noise reference signal. This noise reference signal 

20 is correlated to the undesired signal portion of each of the 
first and second measured signals. 

The noise reference signal is then used to remove the 
undesired portion of each of the first and second measured 
signals via an adaptive noise canceler. An adaptive noise 

25 canceler can be described by analogy to a dynamic multiple 
notch filter which dynamically changes its transfer function 
in response to the noise reference signal and the measured 
signals to remove frequencies from the measured signals that 
are also present in the noise reference signal. Thus, a 

30 typical adaptive noise canceler receives the signal from 
which it is desired to remove noise and a noise reference 
signal. The output of the adaptive noise canceler is a good 
approximation to the desired signal with the noise removed. 
Physiological monitors can often advantageously employ 

35 signal processors of the present invention. Often in 

physiological measurements a first signal comprising a first 
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desired p^^ion and a first undesir^^^>ortion and a second 
signal comprising a second desired portion and a second** 
undesired portion are acquired. The signals may be acquired 
by propagating energy through a patient 1 s body (or a material 
5 which is derived from the body, such as breath, blood, or 

tissue, for example) and measuring an attenuated signal after 
transmission or reflection* Alternatively, the signal may be 
acquired by measuring energy generated by a patient's body, 
such as in electrocardiography. The signals are processed 
10 via the signal processor of the present invention to acquire 
a noise reference signal which is input to an adaptive noise 
canceler. 

One physiological monitoring apparatus which can 
advantageously incorporate the features of the present 

15 invention is a monitoring system which determines a signal 
which is representative of the arterial pulse, called a 
plethysmography wave. This signal can be used in blood 
pressure calculations, blood gas saturation measurements, 
etc. A specific example of such a use is in pulse oximetry 

20 which determines the saturation of oxygen in the blood. In 
this configuration, the desired portion of the signal is the 
arterial blood contribution to attenuation of energy as it 
passes through a portion of the body where blood flows close 
to the skin. The pumping of the heart causes blood flow to 

25 increase and decrease in the arteries in a periodic fashion, 
causing periodic attenuation wherein the periodic waveform is 
the plethysmography waveform representative of the pulse. 

A physiological monitor particularly adapted to pulse 
oximetry oxygen saturation measurement comprises two light 

30 emitting diodes (LED's) which emit light at different 
wavelengths to produce first and second signals. A detector 
registers the attenuation of the two different energy signals 
after each passes through an absorptive media, for example a 
digit such as a finger, or an earlobe. The attenuated 

35 signals generally comprise both desired and undesired signal, 
portions. A static filtering system, such as a band pass 
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filter, removes a portion of the undesired signal which is 
static, or constant, or outside of a known bandwidth of* 
interest, leaving an erratic or random undesired signal 
portion, often caused by notion and often difficult to 
remove, along with the desired signal portion. 

Next, a processor of the present invention removes the 
desired signal portions from the measured signals yielding a 
noise reference signal which is a combination of the 
remaining undesired signal portions. The noise reference 
signal is correlated to both of the undesired signal 
portions. The noise reference signal and at least one of the 
measured signals are input to an adaptive noise canceler 
which removes the random or erratic portion of the undesired 
signal. This yields a good approximation to the desired 
plethysmographic signal as measured at one of the measured 
signal wavelengths. As is known in the art, quantitative 
measurements of the amount of oxygenated blood in the body 
can be determined from the plethysmographic signal in a 
variety of ways. 

One aspect of the present invention is a signal 
processor comprising a detector for receiving a first signal 
which travels along a first propagation path and a second 
signal which travels along a second propagation path wherein 
a portion of the first and second propagation paths are 
located in a propagation medium. The first signal has a 
first desired signal portion and a first undesired signal 
portion and the second signal has a second desired signal 
portion and a second undesired signal portion. The first and 
second undesired signal portions are a result of a 
perturbation of the propagation medium. This aspect of the 
invention additionally comprises a reference processor having 
an input for receiving the first and second signals. The 
processor is adapted to combine the first and second signals 
to generate a reference signal having a primary component 
which is a function of the first and said second undesired 
signal portions. 



The iftve described aspect of the present invention may 
further comprise an adaptive signal processor for receiving — 
the reference signal and the first signal and for deriving 
therefrom an output signal having a primary component which 
5 is a function of the first desired signal portion of the 

first signal. Alternatively, the above described aspect of 
the present invention my further comprise an adaptive signal 
processor for receiving the reference signal and the second 
signal and for deriving therefrom an output signal having a 
10 primary component which is a function of the second desired 
signal portion of the second signal. The adaptive signal 
processor may comprise an adaptive noise canceler. The 
adaptive noise canceler may be comprise a joint process 
estimator having a least-squares-lattice predictor and a 
15 regression filter. 

The detector in the aspect of the signal processor of 
the present invention described above may further comprise a 
sensor for sensing a physiological function. The sensor may 
comprise a light sensitive device. Additionally, the present 
20 invention may further comprising a pulse oximeter for 
measuring oxygen saturation in a living organism. 

Another aspect of the present invention is a 
physiological monitoring apparatus comprising a detector for 
receiving a first physiological measurement signal which 
25 travels along a first propagation path and a second 
physiological measurement signal which travels along a second 
propagation path. A portion of the first and second 
propagation paths is located in a propagation medium. The 
first signal has a first desired signal portion and a first 
30 undesired signal portion and the second signal has a second 
desired signal portion and a second undesired signal portion. 
The physiological monitoring apparatus further comprises a 
reference processor having an input for receiving the first 
and second signals. The processor is adapted to combine the 
35 first and second signals to generate a reference signal 
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having a primary component which is a function of the first 
and the second undesired signal portions. 

The physiological monitoring apparatus may further 
comprise an adaptive signal processor for receiving the 
5 reference signal and the first signal and for deriving 

therefrom an output signal having a primary component which 
is a function of the first desired signal portion of the 
first signal. Alternatively, the physiological monitoring 
apparatus may further comprise an adaptive signal processor 
10 for receiving the reference signal and the second signal and 
for deriving therefrom an output signal having a primary 
component which is a function of the second desired signal 
portion of the second signal. The physiological monitoring 
apparatus may further comprise a pulse oximeter. 
15 A further aspect of the present invention is an 

apparatus for measuring a blood constituent comprising an 
energy source for directing a plurality of predetermined 

' % i wavelengths of electromagnetic energy upon a specimen and a 

detector for receiving the plurality of predetermined 

iz 20 wavelengths of electromagnetic energy from the specimen. The 
detector produces electrical signals corresponding to the 

7' predetermined wavelengths in response to the electromagnetic 

P energy. At least two of the electrical signals each has a 

desired signal portion and an undesired signal portion. 

=ji 25 Additionally, the apparatus comprises a reference processor 
having an input for receiving the electrical signals. The 
processor is configured to combine said electrical signals to 
generate a reference signal having a primary component which 
is derived from the undesired signal portions. 
30 This aspect of the present invention may further 

comprise an adaptive signal processor for receiving the 
reference signal and one of the two electrical signals and 
for deriving therefrom an output signal having a primary 
component which is a function of the desired signal portion 
35 of the electrical signal. This may be accomplished by use of 

an adaptive noise canceler in the adaptive signal processor 
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which may employ a joint process estimator having a least- 
squares-lattice predictor and a regression filter. 

Yet another aspect of the present invention is a blood 
gas monitor for non-invasively measuring a blood constituent 
5 in a body comprising a light source for directing at least 

two predetermined wavelengths of light upon a body and a 
detector for receiving the light from the body. - The 
detector , in response to the light from the body, produces at 
^east two electrical signals corresponding to the at least 

10 two predetermined wavelengths of light. The at least two 
electrical signals each has a desired signal portion and an 
undesired signal portion. The blood oximeter further 
comprises a reference processor having an input for receiving 
the at least two electrical signals. The processor i6 

15 adapted to combine the at least two electrical signals to 
generate a reference signal with a primary component which is 
derived from the undesired signal portions. The blood 
oximeter may further comprise an adaptive signal processor 
for receiving the reference signal and the two electrical 

20 signals and for deriving therefrom at least two output 
signals which are substantially equal, respectively, to the 
desired signal portions of the electrical signals. 

The present invention also includes a method of 
determining a noise reference signal from a first signal 

25 comprising a first desired signal portion and a first noise 
portion and a second signal comprising a second desired 
signal portion and a second noise portion. The method 
comprises the steps of selecting a signal coefficient which 
is proportional to a ratio of predetermined attributes of the 

30 first desired signal portion and predetermined attributes of 

the second desired signal portion. The first signal and the 
second signal coefficient are input into a signal multiplier 
wherein the first signal is multiplied by the signal 
coefficient thereby generating a first intermediate signal. 

35 The second signal and the first intermediate signal are input 

into a signal subtractor wherein the first intermediate 
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signal is subtracted from the second signal. This generates 
a noise reference signal having a primary component which is" 
derived from the first and second noise signal portions. The 
first and second signals in this method may be derived from 
5 light energy transmitted through an absorbing medium. 

The present invention further embodies a physiological 
monitoring apparatus comprising means for acquiring a first 
signal comprising a first desired signal portion and a first 
undesired signal portion and a second signal comprising a 

10 second desired signal portion and a second undesired signal 
portion. The physiological monitoring apparatus of the 
present invention also comprises means for determining from 
the first and second signals a noise reference signal. 
Additionally, the monitoring apparatus comprises an adaptive 

15 noise canceler having a noise reference input for receiving 
the noise reference signal and a signal input for receiving 
the first signal wherein the adaptive noise canceler, in real 
or near real time, generates an output signal which 
approximates the first desired signal portion. The adaptive 

20 noise canceler may further comprise a joint process 
estimator* 

A further aspect of the present invention is an 
apparatus for processing an amplitude modulated signal having 
a signal amplitude complicating feature, the apparatus 

25 comprising an energy source for directing electromagnetic 
energy upon a specimen. Additionally, the apparatus 
comprises a detector for acquiring a first amplitude 
modulated signal and a second amplitude modulated signal. 
Each of the first and second signals has a component 

30 containing information about the attenuation of 
electromagnetic energy by the specimen and a signal amplitude 
complicating feature* The apparatus includes a reference 
processor for receiving the first and second amplitude 
modulated signals and deriving therefrom a noise reference 

35 signal which is correlated with the signal amplitude 
complicating feature. Further, the apparatus incorporates an 
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adaptive noise canceler having a signal input for receiving 
the first amplitude modulated signal, a noise reference input" 
for receiving the noise reference signal, wherein the 
adaptive noise canceler produces an output signal having a 
primary component which is derived from the component 
containing information about the attenuation of 
electromagnetic energy by the specimen. 

Still another aspect of the present invention is an 
apparatus for extracting a plethysmography waveform from an 
amplitude modulated signal having a signal amplitude 
complicating feature, the apparatus comprising a light source 
for transmitting light into an organism and a detector for 
monitoring light from the organism. The detector produces a 
first light attenuation signal and a second light attenuation 
signal, wherein each of the first and second light 
attenuation signals has a component which is representative 
of a plethysmography waveform and a component which is 
representative of the signal amplitude complicating feature. 
The apparatus also includes a reference processor for 
receiving the first and second light attenuation signals and 
deriving therefrom a noise reference signal. The noise 
reference signal and the signal amplitude complicating 
feature each has a frequency spectrum. The frequency 
spectrum of the noise reference signal is correlated with the 
frequency spectrum of the signal amplitude complicating 
feature. Additionally incorporated into this embodiment of 
the present invention is an adaptive noise canceler having a 
signal input for receiving the first attenuation signal and 
a noise reference input for receiving the noise reference 
signal. The adaptive noise canceler produces an output 
signal having a primary component which is derived from the 
component which is representative of a plethysmography 
waveform . 

The present invention also comprises a method of 
removing a motion artifact signal from a signal derived from 
a physiological measurement wherein a first signal having a 
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physiolo^Tcal measurement component and a motion artifact 
component and a second signal having a physiological- 
measurement component and a motion artifact component are 
acquired. From the first and second signals a motion 
artifact noise reference signal which is a primary function 
of the first and second signals motion artifact components is 
deriyed. This method of removing a motion artifact signal 
from a signal derived from a physiological measurement may 
also comprise the step of inputting the motion artifact noise 
reference signal into an adaptive noise canceler to produce 
an output signal which is a primary function of the first 
signal physiological measurement component. 

Brief Description of the Drawings 
Figure 1 illustrates an ideal plethysmographic waveform. 
Figure 2 schematically illustrates the cross-sectional 
structure of a typical finger. 

Figure 3 illustrates a plethysmographic waveform which 
includes a motion-induced undesired erratic signal portion. 

Figure 4 illustrates a schematic diagram of a 
physiological monitor incorporating a processor of the 
present invention and an adaptive noise canceler. 

Figure 4a illustrates the transfer function of a 
multiple notch filter. 

Figure 5 illustrates an example of an adaptive noise 
canceler which could be employed in a physiological monitor 
which also incorporates the processor of the present 
invention. 

Figure 6a illustrates a schematic absorbing material 
comprising N constituents within an absorbing material* 

Figure 6b illustrates another schematic absorbing 
material comprising N constituents within an absorbing 
material. 

Figure 7 is a schematic model of a joint process 
estimator comprising a least-squares lattice predictor and a 
regression filter. 
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Fi^R 8 is a flowchart representing a subroutine 
capable of implementing a joint process estimator as modeled— 
in Figure 7. 

Figure 9 is a schematic model of a joint process 
5 estimator with a least-squares lattice predictor and two 
regression filters. 

Figure 10 is an example of a physiological monitor 
incorporating a processor of the present invention and an 
adaptive noise canceler within a microprocessor. This 
10 physiological monitor is specifically designed to measure a 
plethysmography waveform and perform pulse oximetry 
measurements . 

Figure 11 is a graph of oxygenated and deoxygenated 
absorption coefficients vs. wavelength. 
15 Figure 12 is a graph of the ratio of the absorption 

coefficients of deoxygenated hemoglobin divided by oxygenated 
hemoglobin vs. wavelength. 

Figure 13 is an expanded view of a portion of Figure 11 
marked by a circle labelled 13. 
20 Figure 14 illustrates a signal measured at a first red 

wavelength Aa=Aredl=650nm for use in a processor of the 
present invention employing the ratiometric method for 
determining the noise reference signal n» (t) and for use in 
a joint process estimator. The measured signal comprises a 
25 desired portion Y A .(t) and an undesired portion n A .(t). 

Figure 15 illustrates a signal measured at a second red 
wavelength Ab=Ared2=685nm for use in a processor of the 
present invention employing the ratiometric method for 
determining the noise reference signal n • (t) . The measured 
30 signal comprises a desired portion Y Ab (t) and an undesired 
portion n Ab (t) . 

Figure 16 illustrates a signal measured at an infrared 
wavelength Ac=AlR=940nm for use in a joint process estimator. 
The measured signal comprises a desired portion Y AC (t) and an 
35 undesired portion ni c (t). 
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Fig3?£ 17 illustrates the noisereference signal n' (t) 
determined by a processor of the present invention using the- 
ratiometrlc method. 

Figure 18 illustrates a good approximation Y\.(t) to the 
desired portion Y AB (t) of the signal S A .(t) measured at 
Aa=Aredl=650nm estimated with a noise reference signal n» (t) 
determined by the ratiometric method. 

Figure 19 illustrates a good approximation Y' AC (t) to the 
desired portion Y Ac (t) of the signal s AC (t) measured at 
Ac=AIR=940nm estimated with a noise reference signal n' (t) 
determined by the ratiometric method. 

Figure 20 illustrates a signal Measured at a red 
wavelength Aa=Ared=660nm for use in a processor of the 
present invention employing the constant saturation method 
for determining the noise reference signal n» (t) and for use 
in a joint process estimator. The measured signal comprises 
a desired portion Y A .(t) and an undesired portion n Aa (t). 

Figure 21 illustrates a signal measured at an infrared 
wavelength Ab=AIR=940nm for use in a processor of the present 
invention employing the constant saturation method for 
determining the noise reference signal n» (t) and for use in 
a joint process estimator. The measured signal comprises a 
desired portion Yj^t) and an undesired portion n Ab (t) . 

Figure 22 illustrates the noise reference signal n'(t) 
determined by a processor of the present invention using the 
constant saturation method. 

Figure 23 illustrates a good approximation Y' A ,(t) to the 
desired portion Y,.(t) C f the signal S Al (t) measured at 
Aa=Ared=660nm estimated with a noise reference signal n'(t) 
determined by the constant saturation method. 

Figure 24 illustrates a good approximation Y'^t) to the 
desired portion Y^t) of the signal S^t) measured at 
Ab=AIR=940nm estimated with a noise reference signal n'(t) 
determined by the constant saturation method. 
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Detai T»H n escripfion of the Invent* op 
The present invention is a processor which determines a~" 
noise reference signal n' (t) for use in an adaptive noise 
canceler. An adaptive noise canceler estimates a good 
approximation Y • (t) to a desired signal Y(t) from a composite 
signal S(t) - Y(t) + n(t) which, in addition to the desired 
portion Y(t) comprises an undesired portion n(t) . The 
undesired portion n(t) may contain one or more of a constant 
portion, a predictable portion, an erratic portion, a random 
portion, etc. The approximation to the desired signal Y* (t) 
is derived by removing as many of the undesired portions n(t) 
from the composite signal S(t) as possible. The constant 
portion and predictable portion are easily removed with 
traditional filtering techniques, such as simple subtraction, 
low pass, band pass, and high pass filtering. The erratic 
portion is more difficult to remove due to its unpredictable 
nature. If something is known about the erratic signal, even 
statistically, it could be removed from the measured signal 
via traditional filtering techniques. However, it is often 
the case that no information is known about the erratic 
portion of the noise. In this case, traditional filtering 
techniques are usually insufficient. Often no information 
about the erratic portion of the measured signal is known. 
Thus, an adaptive noise canceler is utilized in the present 
invention to remove the erratic portion. 

Generally, an adaptive noise canceler has two signal 
inputs and one output. One of the inputs is the noise 
reference signal n«(t) which is correlated to the erratic 
undesired signal portions n(t) present in the composite 
signal S (t) . The other input is for the composite signal 
S(t). Ideally, the output of the adaptive noise canceler 
Y'(t) corresponds to the desired signal portion Y(t) only. 
Often, the most difficult task in the application of adaptive 
noise cancelers is determining the noise reference signal 
n'(t) which is correlated to the erratic undesired portion 
n(t) of the measured signal S(t) since, as discussed above, 
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unpredicfSble signal portions are usually quite difficult to 
isolate from the measured signal S(t). In the signal- 
processor of the present invention, a noise reference signal 
n" (t) is determined from two composite signals measured 
simultaneously, or nearly simultaneously, at two different 
wavelengths, Aa and Ab. The signal processor of the present 
invention can be advantageously used in a monitoring device, 
such a monitor being well suited for physiological 
monitoring. 

A block diagram of a generic monitor incorporating a 
signal processor, or reference processor, according to the 
present invention and an adaptive noise canceler is shown in 
Figure 4. Two measured signals, S A .(t) and S^ft) , are 
acquired by a detector 20. One skilled in the art will 
realize that for some physiological measurements, more than 
one detector may be advantageous. Each signal is conditioned 
by a signal conditioner 22a and 22b. Conditioning includes, 
but is not limited to, such procedures as filtering the 
signals to remove constant portions and amplifying the 
signals for ease of manipulation. The signals are then 
converted to digital data by an analog-to-digital converter 
24a and *24b. The first measured signal S AB (t) comprises a 
first desired signal portion, labelled herein Y Aa (t), and a 
first undesired signal portion, labelled herein n Aa (t). The 
second measured signal S&it) is at least partially correlated 
to the first measured signal S Aa (t) and comprises a second 
desired signal portion, labelled herein Y^t), and a second 
undesired signal portion, labelled herein n^Ct) . Typically 
the first and second undesired signal portions, n Aa (t) and 
r*ib(t) r are uncorrelated and/or erratic with respect to the 
desired signal portions Y Aa (t) and Y^t) . The undesired 
signal portions n Aa (t) and n^ft) are often caused by motion 
of a patient. The signals S Aa (t) and S&it) are input to a 
reference processor 26. The reference processor 26 
multiplies the second measured signal S Ab (t) by a factor w and 
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then subtracts the second measured signal S^ft) from the 
first measured signal S Aa (t) . The factor » is determined to 
cause the desired signal portions Y A ,(t) and Y Ab (t) to cancel 
when the two signals S A .(t) and S Ab (t) are subtracted. Thus, 
the output of the reference processor 26 is a noise reference 
signal n« (t) = n Aa (t) - wn^t) which is correlated to both of 
the erratic undesired signal portions n Aa (t) and n^ft) . The 
noise reference signal n'(t) is input, along with one of the 
measured signals Sj^t) , to an adaptive noise canceler 27 
which uses the noise reference signal n' (t) to remove the 
undesired signal portion n Aa (t) or n^t) from the measured 
signal S Aa (t). It will be understood that S^ft) could have 
been input to the adaptive noise canceler 27 along with the 
noise reference signal n • (t) instead of (t) . The output 
of the adaptive noise canceler 27 is a good approximation 
v, Aa(t) to the desired signal Y Aa (t) . The approximation 
v, A«(t) is displayed on the display 28. 

An adaptive noise canceler 30, an example of which is 
shown in block diagram in Figure 5, is employed to remove the 
erratic, undesired signal portions n^t) and n^t) from the 
measured signals S A .(t) and S^t) . The adaptive noise 
canceler 30 in Figure 5 has as one input a sample of the 
noise reference signal n'(t) which is correlated to the 
undesired signal portions n Aa (t) and n^ft) . The noise 
reference signal n' (t) is determined from the two measured 
signals S Aa (t) and S^t) by the processor 26 of the present 
invention as described herein. A second input to the 
adaptive noise canceler is a sample of either the first or 
second measured signal S A .(t) - Y A .(t) + n A .(t) or S^t) «= 
v Ab(t) + n Ab (t) . 

The adaptive noise canceler 30 functions to remove 
frequencies common to both the noise reference signal n*(t) 
and the measured signal S A .(t) or S^t) . since the noise 
reference signal n' (t) is correlated to the erratic undesired 
signal portions n Aa (t) and n Ab (t) , the noise reference signal 
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n'(t) is efft, ; ei-.atic. The adaptive noise nancelsr acts in 
a manner which may be analogized to a dynamic multiple notch 
filter based on the spectral distribution of the noise 
reference signal n'(t). 

Referring to Figure 4a, the transfer function of a 
multiple notch filter is shown. The notches, or dips in the 
amplitude of the transfer function, indicate frequencies 
which are attenuated or removed when a composite measured 
signal passes through the notch filter. The output of the 
notch filter is the composite signal having frequencies at 
which a notch was present removed . In the analogy to an 
adaptive noise canceler, the frequencies at which notches are 
present change continuously based upon the inputs to the 
adaptive noise canceler. 

The adaptive noise canceler 30 shown in Figure 5 
produces an output signal, labelled herein Y^.Ct) or Y f Ab (t), 
which is fed back to an internal processor 32 within the 
adaptive noise canceler 30. The internal processor 32 
automatically adjusts its own transfer function according to 
a predetermined algorithm such that the output of the 
internal processor 32, labelled b(t), closely resembles the 
undesired signal portion n Aa (t) or n^Ct) . The output b(t) of 
the internal processor 32 is subtracted from the measured 
signal, S Aa (t) or S Ab (t), yielding a signal Y» Aa (t) *= S Aa (t) + 
n Aa (t) - b Aa (t) or Y'xbft) « S^t) + nAb (t) - b^t) . The 
internal processor optimizes Y' Aa (t) or Y'^Ct) such that 
* f Aa(t) or Y'xbCt) is approximately equal to the desired 
signal Y Xa (t) or Y^t) , respectively. 

One algorithm which may be used for the adjustment of 
the transfer function of the internal processor 32 is a 
least-squares algorithm, as described in Chapter 6 and 
Chapter 12 of the book Adaptive Signal Processing by Bernard 
Widrow and Samuel Stearns, published by Prentice Hall, 
copyright 1985. This entire book, including Chapters 6 and 
12, is hereby incorporated herein by reference. 



-19- 



Adaptive processors :>0 have been successfully applied to 
a number of problems including antenna sidelobe canceling," 
pattern recognition, the elimination of periodic interference 
in general, and the elimination of echoes on long distance 
telephone transmission lines. However, considerable 

ingenuity is often required to find a suitable noise 
reference signal n'(t) for a _ given application since the 
random or erratic portions n Aa (t) or n^t) cannot easily be 
separated from the measured signal S Aa (t) or S^t) . If the 
actual undesired signal portion n^ft) or n^t) were a priori 
available, techniques such as adaptive noise canceling would 
not be necessary. The unique determination of a suitable 
noise reference signal n'(t) from measurements taken by a 
monitor incorporating a reference processor of the present 
invention is one aspect of the present invention. 

GENERALIZED DETERMINATIO N OF NOISE REFERENCE fiJftpA J- 
An explanation which describes how the noise reference 
signal n' (t) may be determined as follows. A first signal is 
measured at, for example, a wavelength Aa, by a detector 
yielding a signal S Aa (t): 

S A .(t) - Y Aa (t) + n Aa (t); (i) 
where Y Aa (t) is the desired signal and n Aa (t) is the noise 
component . 

A similar measurement is taken simultaneously, or nearly 
simultaneously, at a different wavelength, Ab, yielding: 

Sxb(t) - Y^t) + n Ab . (2 ) 
Note that as long as the measurements, S Aa (t) and S^Ct) , are 
taken substantially simultaneously, the undesired signal 
components, n Aa (t) and n^ft) , will be correlated because any 
random or erratic functions will affect each measurement in 
nearly the same fashion. 

To obtain the noise reference signal n' (t) , the measured 
signals S^t) and S^Ct) are transformed to eliminate the 
desired signal components. One way of doing this is to find 
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a propoz^^nality constant, w,, betveen the desired signals 
Y*i(t) and Y^ft) such that: 

Y>.(t) - v^t) . (3 ) 
This proportionality relationship can be satisfied in many 
5 measurements, including but not limited to absorption 
measurements and physiological measurements. Additionally f 
in most measurements, the proportionality constant can be 
determined such that: 

n Aa (t) h t^njtCt) . (4) 
10 Multiplying equation (2) by w 1 and then subtracting 

equation (2) from equation (1) results in a single equation 
wherein the desired signal terms Y Xm (t) and S^t) cancel, 
leaving: 

n f (t) - S Aa (t) - ti,S»(t) = n Aa (t) - vWt); (5) 
15 a non-zero signal which is correlated to each undesired 
signal portion n A8 (t) and n^t) and can be used as the noise 
reference signal n^t) in an adaptive noise canceler. 

EXAMPLE OF DETERMINATION O F NOISE REFERENCE SIGNAL 
IN AN ABSORPTIVE SYSTEM 

20 Adaptive noise canceling is particularly useful in a 

large number of measurements generally described as 
absorption measurements. An example of an absorption type 
monitor which can advantageously employ adaptive noise 
canceling based upon a noise reference signal n 1 (t) 

25 determined by a processor of the present invention is one 
which determines the concentration of an energy absorbing 
constituent within an absorbing material when the material is 
subject to perturbation. .Such perturbations can be caused by 
forces about which information is desired, or alternatively, 

30 by random or erratic forces such as a mechanical force on the 
material. Random or erratic interference, such as motion, 
generates undesired noise components in the measured signal. 
These undesired components can be removed by the adaptive 
noise canceler if a suitable noise reference signal n 1 (t) is 

35 known . 
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A s^^matic N constituent abso^^ng matarial comprising 
a container 42 having N different absorbing constituents,* 
labelled A,, A^ A3, . A^, is 6hown schematically in Figure 
6a. The constituents A 1 through A^ in Figure 6a are arranged 
5 in a generally orderly, layered fashion within the container 

42. An example of a particular type of absorptive system is 
one in which light energy passes through the container 42 and 
is absorbed according to the generalized Beer-Lambert Law of 
light absorption. For light of wavelength Aa, this 
10 attenuation may be approximated by: 

(6) 

Initially transforming the signal by taking the natural log 
of both sides and manipulating terms, the signal is 
transformed such that the signal components are combined by 
15 addition rather than multiplication, i.e.: 

S A .=ln (I Q /I) - S^^.c^ 

(7) 

where I 0 is the incident light energy intensity; I is the 
transmitted light energy intensity; c c xm is the absorption 
coefficient of the i th constituent at the wavelength Aa; x f (t) 

20 is the optical path length of i th layer, i.e., the thickness 
of material of the i th layer through which optical energy 
passes; and c t (t) is the concentration of the i th constituent 
in the volume associated with the thickness x, (t) . The 
absorption coefficients c, through c r are known values which 

25 are constant at each wavelength. Most concentrations c^t) 
through c M (t) are typically unknown, as are most of the 
optical path lengths x, (t) of each layer. The total optical 
path length is the sum of each of the individual optical path 
lengths x f (t) of each layer. 

30 When the material is not subject to any forces which 

cause perturbation in the thicknesses of the layers, the 
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optical path length of each layer, x f (t), is generally 
constant. This results in generally constant attenuation of 5 
the optical energy and thus, a generally constant offset in 
the measured signal. Typically, this portion of the signal 
5 is of little interest since knowledge about a force which 

perturbs the material is usually desired. Any signal portion 
outside of a known bandwidth of interest, including the 
constant undesired signal portion resulting from the 
generally constant absorption of the constituents when not 

10 subject to perturbation, should be removed* This is easily 
accomplished by traditional band pass filtering techniques. 
However, when the material is subject to forces, each layer 
of constituents may be affected by the perturbation 
differently than each other layer. Some perturbations of the 

15 optical path lengths of each layer x f (t) may result in 
excursions in the measured signal which represent desired 
information. Other perturbations of the optical path length 
of each layer x ( (t) cause undesired excursions which mask 
desired information in the measured signal. Undesired signal 

20 components associated with undesired excursions must also be 
removed to obtain desired information from the measured 
signal. 

The adaptive noise canceler removes from the composite 
signal, measured after being transmitted through or reflected 

25 from the absorbing material, the undesired signal components 
caused by forces which perturb the material differently from 
the forces which perturbed the material to cause the desired 
signal component. For the purposes of illustration, it will 
be assumed that the portion of the measured signal which is 

30 deemed the desired signal Y At (t) is the attenuation term 
£505X5 (t) associated with a constituent of interest, namely A^ 
and that the layer of constituent Aj is affected by 
perturbations differently than each of the layers of other 
constituents A t through A 4 and A 6 through A,,. An example of 

35 such a situation is when layer A^ is subject to forces about 
which information is desired and, additionally, the entire 
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material is subject to forces which affect each of the 
layers. In this case, since the total force affecting the" 
layer of constituents is different than the total forces 
affecting each of the other layers and information is desired 
about the forces and resultant perturbation of the layer of 
constituents Aj, attenuation terms due to constituents A, 
through A 4 and A 6 through A^ make up the undesired signal 
n A .(t). Even if the additional forces which affect the entire 
material cause the same perturbation in each layer, including 
the layer of Aj, the total forces on the layer of constituent 
Aj cause it to have different total perturbation than each of 
the other layers of constituents A 1 through A 4 and A 6 through 
V 

It is often the case that the total perturbation 
affecting the layers associated with the undesired signal 
components is caused by random or erratic forces. This 
causes the thickness of layers to change erratically and the 
optical path length of each layer, x, (t) , to change 
erratically, thereby producing a random or erratic undesired 
signal component n^(t). However, regardless of whether or 
not the undesired signal portion n A ,(t) is erratic, the 
undesired signal component n Aa (t) can be removed via an 
adaptive noise canceler having as one input a noise reference 
signal n^t) determined by a processor of the present 
invention as long as the perturbation on layers other than 
the layer of constituent Aj is different than the 
perturbation on the layer of constituent Aj. The adaptive 
noise canceler yields a good approximation to the desired 
signal Y' Aa (t). From this approximation, the concentration 
of the constituent of interest, CjCt), can often be 
determined since in some physiological measurements, the 
thickness of the desired signal component, Xj(t) in this 
example, is known or can be determined. 

The adaptive noise canceler utilizes a sample of a noise 
reference signal n'(t) determined from two substantially 
simultaneously measured signals S Aa (t) and S Ab (t) . S Aa (t) is 
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determined as above in equation (7). S xb (t) is determined 
similarly at a different wavelength xb. To find the noise 
reference signal n'(t), attenuated transmitted energy is 
measured at the two different wavelengths Aa and xb and 
transformed via logarithmic conversion. The signals S x ,(t) 
and S^bCt) can then be written (logarithm converted) as: 



S u ( t) - e StU c sXs ( t) ♦ I ^c-jjcj ( t) ♦E£«e jt<Xa c Jc x Jt ( t) ] 



«s.Aa<V^(t)+ n u (t) 



5^( t) = e 5tXb c sXs ( t) ♦ [Si-ie^c^ ( t) *2£«e Jt#U) c it x Jt ( t) ] 



(8) 



(9) 
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(11) 

A further transformation of the signals is the 
proportionality relationship defining ug, similarly to 
equation (3) , which allows determination of a noise reference 
signal n' (t) , is: 

c s.xm = «z e 5.Ab* ^ere (12) 
n x* * **n*b- (13) 
It is often the case that both equations (12) and (13) can be 
simultaneously satisfied. Multiplying equation (11) by 
and subtracting the result from equation (9) yields a non- 
20 zero noise reference signal which is a linear sum of 
undesired signal components: 

n»(t) = S A .(t) - WjS^ft) m n A ,(t) - ^ Ab (t). (14) 

- 2 '-i e i. x* c i x i < c > ♦SJU«*. u^iXk ( t) 
-El. 1 co 2 e i t j^c^ ( t) a 2 e kt ( c) 

• (15) 
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(16) 

A sample of this noise reference signal n'(t), and a 
sample of either measured signal S A ,(t) or S ib (t) , are input 
to an adaptive noise canceler, one model of which is shown in 
Figure 5 and a preferred .model of which is discussed herein 
under the heading PREFERRED ADAPTIVE NOISE CANCELER USING A 
JOINT PROCESS ESTIMATOR IMPLEMENTATION. The adaptive noise 
canceler removes the undesired portion of the measured signal 
n*«(t) or nib(t) , yielding a good approximation to the desired 
portion of signal Y',.(t) « « 5<u <ycs(t) or Y'^t) - c^cys, (t) . 
The concentration Cj(t) may then be determined from the 
approximation to the desired signal Y' A ,(t) or Y'^ft) 
according to: 

c 5 (t) « V'x.WAs^ft) . Y' Ab (t)/e 5fXbX5 (t). (17) 
As discussed previously, the absorption coefficients are 
constant at each wavelength Aa and Ab and the thickness of 
the desired signal component, Xj(t) in this example, is often 
known or can be determined as a function of time, thereby 
allowing calculation of the concentration c^(t) of 
constituent Aj. 

DETERMINATION OF CONCENTRATION OR SATURATION 
IN A VOLUME CONTAINING MORE THAN ONE CONSTITUENT 
Referring to Figure 6b, another material having N 
different constituents arranged in layers is shown. In this 
material, two constituents Aj and Aj are found within one 
layer having thickness Xj 6 (t) - x^t) + x 6 (t) located 
generally randomly within the layer. This is analogous to 
combining the layers of constituents A, and A 6 in Figure 6a. 
A combination of layers, such as the combination of layers of 
constituents Aj and A^, is feasible when the two layers are 
under the same total forces which result in the same 
perturbation of the optical path lengths Xj(t) and x 6 (t) of 
the layers. 




Often it is desirable to find the concentration or the 
saturation, i.e., a percent concentration, of one constituent- 
within a given thickness which contains more than one 
constituent and is subject to unique forces. A determination 
5 of the concentration or the saturation of a constituent 

within a given volume may be made with any number of 
constituents in the volume subject to the same total forces 
and therefore under the same perturbation. To determine the 
saturation of one constituent in a volume comprising many 
10 constituents, as many measured signals as there are 

constituents which absorb incident light energy are 
necessary. It will be understood that constituents which do 
not absorb light energy are not consequential in the 
determination of saturation. To determine the concentration, 
15 as many signals as there are constituents which absorb 
incident light energy are necessary as well as information 
about the sum of concentrations. 

It is often the case that a thickness under unique 
motion contains only two constituents. For example, it may 
20 be desirable to know the concentration or saturation of Aj 
within a given volume which contains Aj and A 6 . In this case, 
the desired signals Y A- (t) and Y^t) comprise terms related 
to both Aj and A 6 so that a determination of the concentration 
or saturation of A5 or A 6 in the volume may be made. A 
25 determination of saturation is discussed herein. It will be 
understood that the concentration of Aj in volume containing 
both A5 and A 6 could also be determined if it is known that 
A5+A 6 =l, i.e., that there are no constituents in the volume 
which do not absorb incident light energy at the particular 
30 measurement wavelengths chosen. The measured signals S Ai (t) 
and S Ab (t) can be written (logarithm converted) as: 

S A .(t) = £5^.05X5 6 (t) + e 6#i .c 6 X5 6 (t) + n A .(t) (18) 
- Y A .(t) + n A .(t) ; (19) 
S Ab (t) = cs.^csx^t) + c 6 . Ab c 6 X5 6 (t) + n xb (t) (20) 
35 "* ** b (t) + n Ab (t). (21) 
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Any signal portions outside of a known bandwidth of 
interest, including the constant undesired signal portiorr 
resulting from the generally constant absorption of the 
constituents when not under perturbation, should be removed 
5 to determine an approximation to the desired signal. This is 

easily accomplished by traditional band pass filtering 
techniques. As in the previous example , it is often the case 
that the total perturbation affecting the layers associated 
with the undesired signal components is caused by random or 
10 erratic forces , causing the thickness of each layer, or the 
optical path length of each layer # x f (t), to change 
erratically, producing a random or erratic undesired signal 
component n Aa (t). Regardless of whether or not the undesired 
signal portion n Aa (t) is erratic, the undesired signal 
15 component n Att (t) can be removed via an adaptive noise canceler 
having as one input a noise reference signal n' (t) determined 
by a processor of the present invention as long as the 
perturbation in layers other than the layer of constituents 
Aj and A 6 is different than the perturbation in the layer of 
20 constituents Aj and A 6 « The erratic undesired signal 
components n Aa (t) and n^t) may advantageously be removed 
from equations (18) and (20), or alternatively equations (19) 
and (21) , by an adaptive noise canceler. The adaptive noise 
canceler, again, requires a sample of a noise reference 
25 signal n' (t) . 

DETERMINATION OF NOISE REFERENCE SIGNAL 
FOR SATURATION MEASUREMENT 

Two methods which may be used by a processor of the 
present invention to determine the noise reference signal 

30 n f (t) are a ratidmetric method and a constant saturation 
method. The preferred embodiment of a physiological monitor 
incorporating a processor of the present invention utilizes 
the ratiometric method wherein the two wavelengths Aa and Ab, 
at which the signals S Aa (t) and S Ab (t) are measured, are 

35 specifically chosen such that a relationship between the 

absorption coefficients t s§Xa , c Sxb , c 6xa and c 6xb exists, i.e.: 



£ S.i»/ c 6.i. ~ «S.*t/ e «.»b (22) 

The measured signals s A .(t) and S^t) can be factored" 
and written as: 

S >^ fc ) " c «. i .C< e S. i ./ e 6. i .> c 3*( t > + c 6 x(t)J + n A8 (t) (23 ) 
S ib (t) - ^[(e^c^JCjXCt) + c 6 x(t)] + n Ab (t). (24) 
The wavelengths Aa and Ab, chosen to satisfy, equation (22) 
cause the terms within the square brackets to" be equal,' 
thereby causing the desired signal portions Y' A .(t) and 
. Y'^ft) to be linearly dependent. Then, a proportionality 
constant w rZ which causes the desired signal portions Y' A ,(t) 
and Y'^ft) to be equal and allows determination of a non-zero 
noise reference signal n'(t) is: 

c 6.a. * v rz%xb' wt > e re (2 5) 
«a« »* w r s n Ab» (26) 
It is often the case that both equations (25) and (26) can be 
simultaneously satisfied. Additionally, since absorption 
coefficients of each constituent are constant with respect to 
wavelength, the proportionality constant „ r3 can be easily 
determined. Furthermore, absorption coefficients of other 
constituents A, through A, and A, through A, are generally 
unequal to the absorption coefficients of Aj and A 6 . Thus, 
the undesired noise components n Att and are generally not 
made linearly dependent by the relationships of equations 
(22) and (25) . 

Multiplying equation (24) by « r3 and subtracting the 
resulting equation from equation (23), a non-zero noise 
reference signal is determined by: 

»'(t) - S A .(t) - t^S^t) = n Ai (t) - « r3nA5 (t). (27) 

An alternative method for determining the noise 
reference signal from the measured signals S A .(t) and S Ab (t) 
using a processor of the present invention is the constant 
saturation approach. In this approach, it is assumed that 
the saturation of % in the volume containing ^ and A 6 
remains relatively constant, i.e.: 

Saturation (^ (t)) = c 5 (t)/[c 5 (t) + c 6 (t)) (28) 



- (1 + [c 6 (t) /C5 (t)])-' ( 
is substantially constant over many samples of the measured- 
signals s xa and S Ab . This assumption is accurate over many 
samples since saturation generally changes relatively slovlv 
5 in physiological systems. 7 

The constant saturation assumption is equivalent to 
assuming that: 

c 5( t )/c 6 (t) - constant (30) 
since the only other term in equation (29) is a constant, 
10 namely the numeral 1. 

Using this assumption/ the proportionality constant 
•*(t) which allows determination of the noise reference 
signal n' (t) is: 

Q,j(t)= g s.u<yfr. t ( t) ♦ e^c^ At) 
€ s.Ju,<Vr 5#< <t) + e 6iil> c 6 ^ 6 (t) 



15 



(31) 

- ^.(tJ/Y^t) (32) 

e S.lb C 5 * e 6.lb C € 



(33) 



20 



(34) 

15 v, A«(t)/Y' Ab (t) - constant; where (35) 
n>.(t) * «- d (t) nAb (t). (36) 
It is often the case that both equations (35) and (36) can be 
simultaneously satisfied to determine the proportionality 
constant ^(t). Additionally, the absorption coefficients 

at each wavelength c r . * , , 

' 5.A8' c 6,xf *5,xb' and « 6 ^ are constant and 

the central assumption of the constant saturation method is 
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that C3(t)/c 6 (t) is constant over many sample periods. Thus 
a new proportionality constant w 0 (t) may be determined every- 
few samples from new approximations to the desired signal as 
output from the adaptive noise canceler. Thus, the 
approximations to the desired signals Y'^t) and Y^ft) 
found by the adaptive noise canceler for a substantially 
immediately preceding set of samples of the measured signals 
S A ,(t) and S Ab (t) are used in a processor of the present 
invention for calculating the proportionality constant, 
w o( t )» for the next set of samples of the measured signals 
S Aa (t) and Sj^t) . 

Multiplying equation (20) by «^(t) and subtracting the 
resulting equation from equation (18) yields a non-zero noise 
reference signal: 

n'(t) = S A .(t) - « IJ (t)S ib (t) = n^t) - w^(t)n Ab (t). (37) 
It will be understood that equation (21) could be multiplied 
by w,s(t) and the resulting equation could be subtracted from 
equation (19) to yield the same noise reference signal n' (t) 
as given in equation (37). 

When using the constant saturation method, it is 
necessary for the patient to remain motionless -for -a short 
period of time such that an accurate initial saturation value 
can be determined by known methods other than adaptive noise 
canceling on which all other calculations will be based. 
With no erratic, motion-induced undesired signal portions, a 
physiological monitor can very quickly produce an initial 
value of the saturation of Aj in the volume containing Aj and 
A 6 . An example of a saturation calculation is given in the 
article "SPECTROPHOTOMETRY DETERMINATION OF OXYGEN 
SATURATION OF BLOOD INDEPENDENT OF THE PRESENT OF INDOCYANINE 
GREEN" by G.A. Mook, et al., wherein determination of oxygen 
saturation in arterial blood is discussed. Another article 
discussing the calculation of oxygen saturation is "PULSE 
OXIMETRY: PHYSICAL PRINCIPLES, TECHNICAL REALIZATION AND 
PRESENT LIMITATIONS" by Michael R. Neuman. Then, with values 



for Y' A8 ^and Y\ b (t) determined, Adaptive noise canceler 
may be utilized with a noise reference signal ^re- 
determined by the constant saturation method. 

PREFERRED ADAPTTVfl p n ISE cANrFT.TSR USING ft 
JQIWT PROCESS ESTIMATOR TM PLEMENTATTOM 
Once the noise reference signal n'(t) is determined by 
the processor of the present invention using either the above 
described ratiometric or constant saturation methods, the 
adaptive noise canceler can be implemented in either hardware 
br software. 

The least mean squares (i*s) implementation of the 
internal processor 32 described above in conjunction with the 
adaptive noise canceler of Figure 5 is relatively easy to 
implement, but lacks the speed of adaptation desirable for 
most physiological monitoring applications of the present 
invention. Thus, a faster approach for adaptive noise 
canceling, called a least-squares lattice joint process 
estimator model, is preferably used. A joint process 
estimator 60 is shown diagrammatical ly in Figure 7 and is 
described in detail in Chapter 9 of Adaptive f<u w T h-~rY by 
Simon Haykin, published by Prentice-Hall, copyright 1986. 
This entire book, including Chapter 9/ is hereby incorporated 
herein by reference. The function of the joint process 
estimator is to remove the undesired signal portions n Aa (t) 
or n Ab (t) from the measured signals S A .(t) or S^(t), yielding 
a signal Y' AB (t) or Y'^t) which is a good approximation to 
the desired signal Y^t) or Y^t) . Thus, the joint process 
estimator estimates the value of the desired signal Y^t) or 
Y^t). The inputs to the joint process estimator 60 are the 
noise reference signal n'(t) and the composite measured 
signal S A .(t) or S^t) . The output is a good approximation 
to the signal s^t) or S*(t) with the noise removed, i.e. a 
good approximation to Y A ,(t) or Y^t) . 

The joint process estimator 60 utilizes, in conjunction, 
a least square lattice predictor 70 and a regression filter 
80. The noise reference signal n» (t) is input to the least 
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square lattice predictor 70 while the measured signal s,.( t ) 
or S^it) is input to the regression filter 80. " F or" 
simplicity in the following description, S Aa (t) will be the 
measured signal from which the desired portion Y ia (t) will be 
estimated by the joint process estimator 60. However, it 
will be noted that S^ft) could equally well be input to' the 
regression filter 80 and the desired portion Y^t) of l^is 
signal could equally well be estimated. 

The joint process estimator 60 removes all frequencies 
that are present in both the noise reference signal n'(t) and 
the measured signal S A .(t). The undesired signal portion 
n Aa (t) usually comprises frequencies unrelated to those of the 
desired signal portion Y A .(t). It is highly improbable that 
the undesired signal portion n A .(t) would be of exactly the 
same spectral content as the desired signal portion Y Aa (t) . 
However, in the unlikely event that the spectral content of 
S A .(t) and n'(t) are similar, this approach will not yield 
accurate results. Functionally, the joint process estimator 
60 compares input signal n»(t), which is correlated to the 
undesired signal portion n Aa (t) , and input signal S Aa (t) and 
removes all frequencies which are identical. Thus, the joint 
process estimator 60 acts as a dynamic multiple notch filter 
to remove those frequencies in the undesired signal component 
n A .(t) as they change erratically with the motion of the 
patient. This yields a signal having substantially the same 
spectral content as the desired signal Y A .(t) . The output of 
the joint process estimator 60 has substantially the same 
spectral content and amplitude as the desired signal Y A .(t) . 
Thus, the output Y' Al (t) of the joint process estimator 60 is 
a very good approximation to the desired signal Y Aa (t). 

The joint process estimator 60 can be divided into 
stages, beginning with a zero-stage and terminating in an up- 
stage, as shown in Figure 7. Each stage, except for the 
zero-stage, is identical to every other stage. The zero- 
stage is an input stage for the joint process estimator 60. 
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The first stage through the m« h -stage work on the signal 
produced in the immediately previous stage, ile., the 
stage, such that a good approximation to the desired signal 
Y'ia(t) is produced as output from the m th -stage. 

The least-squares lattice predictor 70 comprises 
registers 90 and 92, summing elements 100 and 102, and delay 
elements HOr The registers 90 and 92 contain multiplicative 
values of a forward reflection coefficient r f B (t) and a 
backward reflection coefficient r bm (t) which multiply the 
noise reference signal n'(t) and signals derived from the 
noise reference signal n' (t) . Each stage of the least- 
squares lattice predictor outputs a forward prediction error 
f.(t) and a backward prediction error b.(t) . The subscript 
m is indicative of the stage. 

For each set of samples, i.e. one sample of the noise 
reference signal n'(t) derived substantially simultaneously 
with one sample of the measured signal S A .(t), the sample of 
the noise reference signal n'(t) is input to the least- 
squares lattice predictor 70. The zero-stage forward 
prediction error f 0 (t) and the zero-stage backward prediction 
error b 0 (t) are set equal to the noise reference signal 
n» (t) . The backward prediction error b 0 (t) is delayed by one 
sample period by the delay element lio in the first stage of 
the least-squares lattice predictor 70. Thus, the 
immediately previous value of the noise reference signal 
n' (t) is used in calculations involving the first-stage delay 
element lio. The zero-stage forward prediction error is 
added to the negative of the delayed zero-stage backward 
prediction error b 0 (t-l) multiplied by the forward reflection 
coefficient value r f#1 (t) register 90 value, to produce a 
first-stage forward prediction error f , (t) . Additionally, 
the zero-stage forward prediction error f 0 (t) is multiplied 
by the backward reflection coefficient value r b 1{t) register 
92 value and added to the delayed zero-stage backward 
prediction error b 0 (t-l) to produce a first-stage backward 
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prediction error b,(t). in each subsequent stage, m, of the 
least square lattice predictor 70, the previous forward and- 
backward prediction error values, f^,(t) and b^t-l), the 
backward prediction error being delayed by one sample period, 
are used to produce values of the forward and backward 
prediction errors for the present «tage, f B (t) and b B (t) . 

The backward prediction error b.(t) is fed to the 
concurrent stage, m, of the regression filter 80. There it 
is input to a register 96, which contains a multiplicative 
regression coefficient value K m XM (t) . For -example, in the 
zero-stage of the regression filter 80, the rero-stage 
backward prediction error b 0 (t) is multiplied by the zero- 
stage regression coefficient * 0 Aa (t) register 96 value and 
subtracted from the measured value of the signal S Aa (t) at a 
summing element 106 to produce a first stage estimation error 
signal e 1ia (t). The first-stage estimation error signal 
e i,*a( t ) is a first approximation to the desired signal. This 
first-stage estimation error signal e, Xa (t) is input to the 
first-stage of the regression filter 80. The first-stage 
backward prediction error b,(t), multiplied by the first- 
stage regression coefficient * 1<XB (t) register 96 value is 
subtracted from the first-stage estimation error signal 
e i.Aa( t ) to produce the second-stage estimation error e 2xs (t). 
The second-stage estimation error signal e 2 , a (t) is a second, 
somewhat better approximation to the desired signal Y A ,(t). 

The same processes are repeated in the least-squares 
lattice predictor 70 and the regression filter 80 for each 
stage until a good approximation to the desired signal Y» A ,(t) 
= e..*.^) is determined. Each of the signals discussed 
above, including the forward prediction error f.(t) , the 
backward prediction error b m (t) , the estimation error signal 
^.(tK is necessary to calculate the forward reflection 
coefficient r f „(t) , the backward reflection coefficient 
^.■(t) > and the regression coefficient /c - Aa (t) register 90, 
92, and 96 values in each stage, m. In addition to the 



forward prediction error f.(t) , the backward prediction error 
b.(t) , and the estimation error e. x .(t) signals, a number of. 
intermediate variables, not shown in Figure 7 but based on 
the values labelled in Figure 7, are required to calculate 
the forward reflection coefficient r fB (t) , the backward 
reflection coefficient r b .(t) , and the regression coefficient 
V»« (t) re 9ister 90,92, and 96 values. 

Intermediate variables include a weighted sum of the 
forward prediction error squares W m {t) , a weighted sun of the 
backward prediction error squares fi m (t) , a scaler parameter 
A,(t) , a conversion factor 7- (t) , and another scaler parameter 
fim.xmt*-) • The weighted sum of the forward prediction errors 
^.(t) is defined as: 

where x without a wavelength identifier, a or b, is a 
constant multiplicative value unrelated to wavelength and is 
typically less than or equal to one, i.e., x s l. The 
weighted sum of the backward prediction errors fi m (t) is 
defined as: 

t - " 

0 (t) = Z A*** |b (i)| 2 (39) 
m i=1 m 

where, again, X without a wavelength identifier, a or b, is 
a constant multiplicative value unrelated to wavelength and 
is typically less than or equal to one, i.e., Ail. These 
weighted sum intermediate error signals can be manipulated 
such that they are more easily solved for, as described in 
Chapter 9, § 9.3. and defined hereinafter in equations (53) 
and (54). 

DESCRIPTION OF THE JOINT PROCES S ESTIMATOR 
The operation of the joint process estimator 60 is as 
follows. When the joint process estimator 60 is turned on, 
the initial values of intermediate variable and signal 
including the parameter (t) , the weighted sum of the 



forward prediction error signals ^(t), the weighted sum of 
the backward prediction error signals ^(t), the parameter" 
Pm.xal*) * and the zero-stage estimation error e„ A- (t) are 
initialized, some to zero and some to a small positive number 
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Vi(°) - 0; (40) 

= H ' (41) 

/WO) = Si (42) 

Pm.x.(°) " 0? (43) 
^.(t) - S^(t) for t Js 0. (44) 
After initialization, a simultaneous sample of the 
measured signal S A ,(t) and the noise reference signal n«(t) 
are input to the joint process estimator 60, as shown in 
Figure 7. The forward and backward prediction error signals 
f 0 (t) and b 0 (t), and intermediate variables including the 
weighted sums of the forward and backward error signals 7 0 (t) 
and p 0 (t) t and the conversion factor 7 0 (t) are calculated for 
the zero-stage according to: 

M*) - Mt) *= n'(t) (45) 
Mt> - ^o(t) - AM*" 1 ) + I n'(t)l 2 (46) 
7 0 (t-l) = 1 • (47 , 

where, again, x without a wavelength identifier, a or b, is 
a constant multiplicative value unrelated to wavelength. 

Forward reflection coefficient r f .(t) , backward 
reflection coefficient r b .(t), and regression coefficient 
register 90, 92 and 96 values in each stage thereafter 
are set according to the output of the previous stage. The 
forward reflection coefficient r fl (t), backward reflection 
coefficient r b1 (t), and regression coefficient Aa (t) 
register 90, 92 and 96 values in the first stage are thus set 
according to algorithm using values in the zero-stage of the 
joint process estimator 60. In each stage, msl, 

intermediate values and register values including the 
parameter VtW? the forward reflection coefficient r f B (t) 
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register 90 value; the backward reflection coefficient r b- (t) 
register 92 value; the forward and backward error signals" 
f„(t) and b^t) ; the weighted sum of squared forward 
prediction errors * f# .(t) , as manipulated in § 9.3 of the 
Haykin book; the weighted sum of squared backward prediction 
errors 0 bm (t) , as manipulated in § 9.3 of the Haykin book; 
the conversion factor 7- (t) ; the parameter P B<i .(t); the 
regression coefficient K_^(t) register 96 value; and the 
estimation error e^ l Aa (t) value are set according to: 

Vi( fc ) " AVi^" 1 ) + {^(t-lJfVi^JA^Ct-l)} (48) 

T f.ml*) = -(^(^/^(t-l)} (49) 

"-^ViW/'hW) (50) 
f *(V - f^i(t) + r% .(t)b^(t-l) ( 51) 
b.(t) ■= b^t-l) + ^(tj^tt) (52) 
15 V*) - ^i(t) " (I Vl(t)l V^,(t-1)) (53) 

PmW = ^i(t-l) - {I V^.,(t)) (54) 

7.(t-l) = 7»,(t-l) - (I b^ct-Dl V^,(t-D) (55) 
- V., A .(t-l) + {b.(t)e* B Aa (t)/ 7|| (t)) (56) 
= {/>..A.(t)//>.(t)) (57) 
e ~i.*.(t> = e^ x .(t) - *\(t)b m {t) 
where a (*) denotes a complex conjugate. 

These equations cause the error signals f„(t) , b a (t) , 
^■(t) to be squared or to be multiplied by one another, in 
effect squaring the errors, and creating new intermediate 
25 error values, such as Vi( t >- The error signals and the 
intermediate error values are recursively tied together, as 
shown in the above equations (48) through (58). They 
interact to minimize the error signals in the next stage. 
After a good approximation to the desired signal Y' Aa (t) 
30 has been determined by the joint process estimator 60, a next 
set of samples, including a sample of the measured signal 
S Aa (t) and a sample of the noise reference signal n'(t), are 
input to the joint process estimator 60. The re- 
initialization, process does not re-occur, such that the 



(58) 
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forward and backward reflection coefficient r f m (t) and r bi| (t) 
register 90, 92 values and the regression coefficient « a u (t)~ 
register 96 value reflect the multiplicative values required 
to estimate the desired portion Y A .(t) of the sample of S At (t) 
input previously. Thus, information from previous samples is 
used to estimate the desired signal portion of a present set 
of samples in each stage. 

FLOWCHART OF JOINT PROCESS ESTIMATOR 
In a signal processor, such as a physiological monitor, 
incorporating a reference processor of the present invention 
to determine a noise reference signal n f (t) for input to an 
adaptive noise canceler, a joint process estimator 60 type 
adaptive noise canceler is generally implemented via a 
software program having an iterative loop. One iteration of 
the loop is analogous to a single stage of the joint process 
estimator as shown in Figure 7. Thus, if a loop is iterated 
m times, it is equivalent to an m 6tage joint process 
estimator 60. 

A flow chart of a subroutine to estimate the desired 
signal portion Y Aa (t) of a sample of a measured signal, S^t) 
is shown in Figure 8. The flow chart describes how the 
action of a reference processor for determining the noise 
reference signal and the joint process estimator 60 would be 
implemented in software. 

A one-time only initialization is performed when the 
physiological monitor is turned on, as indicated by an 
"INITIALIZE NOISE CANCELER" box 120. The initialization sets 
all registers 90, 92, and 96 and delay element variables 110 
to the values described above in equations (40) through (44). 

Next, a set of simultaneous samples of the measured 
signals S xm (t) and S^Jt) is input to the subroutine 
represented by the flowchart in Figure 8. Then a time update 
of each of the delay element program variables occurs, as 
indicated in a "TIME UPDATE OF [Z* 1 ] ELEMENTS" box 130, 
wherein the value -stored in each of the delay element 
variables 110 is set to the value at the input of the delay 
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element variable 110. Thus, the zero-stage backward 
prediction error b 0 (t) is stored in the first-stage delay- 
element variable, the first-stage backward prediction error 
b,(t) is stored in the second-stage delay element variable, 
and so on. 

Then, using the set of measured signal samples S A ,(t) and 
s Ab(t) , the noise reference signal is calculated according to 
the ratiometric or the constant saturation method described 
above. This is indicated by a "CALCULATE NOISE REFERENCE 
(n'(t)) FOR TWO MEASURED SIGNAL SAMPLES" box 140. The 
ratiometric method is generally preferred since no 
assumptions about constant saturation values need be made. 

A zero-stage order update is performed next as indicated 
in a "ZERO-STAGE UPDATE" box 150. The zero-stage backward 
prediction error b 0 (t), and the zero-stage forward prediction 
error f 0 (t) are set equal to the value of the noise reference 
signal n«(t). Additionally, the weighted sum of the forward 
prediction errors 7 m (t) and the weighted sum of backward 
prediction errors p m (t) are set equal to the value defined in 
equation (46) . 

Next, a loop counter, m, is initialized as indicated in 
a "m=0" box 160. A maximum value of m, defining the total 
number of stages to be used by the subroutine corresponding 
to the flowchart in Figure 8, is also defined. Typically, 
the loop is constructed such that it stops iterating once a 
criterion for convergence upon a best approximation to the 
desired signal has been met by the joint process estimator 
60. Additionally, a maximum number of loop iterations may be 
chosen at which the loop stops iteration. In a preferred 
embodiment of a physiological monitor of the present 
invention, a maximum number of iterations, m=60 to m=80, is 
advantageously chosen. 

Within the loop, the forward and backward reflection 
coefficient r f .(t) and r b B (t) register 90 and 92 values in the 
least-squares lattice filter are calculated first, as 




indicated by the "ORDER UPDATE MTH CELL OF LSL- LATTICE" box 
170 in Figure 8. This requires calculation of intermediate- 
variable and signal values used in determining register 90, 
92, and 96 values in the present stage, the next stage, and 
5 in the regression filter 80 

The calculation of regression filter register 96 value 
K m A.ft) * s Performed next, indicated by the "ORDER UPDATE MTH 
STAGE OF REGRESSION FILTER(S)" box 180. The two order update 
.boxes 170 and 180 are performed in sequence m times, until m 

10 has reached its predetermined maximum (in the preferred 
embodiment, m=60 to m=80) or a solution has been converged 
upon, as indicated by a YES path from a "DONE" decision box 
190. In a computer subroutine, convergence is determined by 
checking if the weighted sums of the forward and backward 

15 prediction errors 7 m (t) and ^(t) are less than a small 
positive number. An output is calculated next, as indicated 
by a "CALCULATE OUTPUT" box 200. The output is a good 
approximation to the desired signal, as determined by the 
reference processor and joint process estimator 60 subroutine 

20 corresponding to the flow chart of Figure 8. This is 
displayed (or used in a calculation in another subroutine) , 
as indicated by a "TO DISPLAY" box 210. 

A new set of samples of the two measured signals S Aa (t) 
and S Ab (t) is input to the processor and joint process 

25 estimator 60 adaptive noise canceler subroutine corresponding 
to the flowchart of Figure 8 and the process reiterates for 
these samples. Note, however, that the initialization 
process does not re-occur. New sets of measured signal 
samples S Aa (t) and S Ab (t) are continuously input to the 

30 reference processor and joint process estimator 60 adaptive 
noise canceler subroutine. The output forms a chain of 
samples which is representative of a continuous wave. This 
waveform is a good approximation to the desired signal 
waveform Y • Aa ( t ) at wavelength Aa. 
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CALCULATION OF SATURATION F ROM AD^PTTVE NOISE CANCELER OUfp rrP 

Physiological monitors typically use the approximation 

of the desired signal Y' A .(t) to calculate another quantity, 
such as the saturation of one constituent in a volume 
5 containing that constituent plus one or more other 
constituents. Generally, such calculations require 

information about a desired signal at two wavelengths. in 
some measurements, this wavelength is Ab, the wavelength used 
in the calculation of the noise reference signal n'(t). For 

10 example, the constant saturation method of determining the 
noise reference signal requires a good approximation of the 
desired signal portions Y Aa (t) and Yj^t) of both measured 
signals S Aa (t) and S^t) . Then, the saturation is determined 
from the approximations to both signals, i.e. Y« Aa (t) and 

15 Y'^ft). 

In other physiological measurements, information about 
a signal at a third wavelength is necessary. For example, to 
find the saturation using the ratiometric method, signals 
S Aa (t) and Sxb(t) are used to find the noise reference signal 

20 n'(t). But as discussed previously, Aa and Ab were chosen to 
satisfy a proportionality relationship like that of equation 
(22). This proportionality relationship forces the two 
desired signal portions Y Aa (t) and Y^t) to be linearly 
dependant. Generally, linearly dependant mathematical 

25 equations cannot be solved for the unknowns. Analogously, 
some desirable information cannot be derived from two 
linearly dependent signals. Thus, to determine the 
saturation using the ratiometric method, a third signal is 
simultaneously measured at wavelength Ac. The wavelength Ac 

30 is chosen such that the desired portion Y Ae (t) of the measured 
signal S Ac (t) is not linearly dependent with the desired 
portions Y Aa (t) and Y^t) of the measured signals S^t) and 
S*b(t) . Since all measurements are taken substantially 
simultaneously, the noise reference signal n'(t) is 

35 correlated to the undesired signal portions n Aa , n Ab , and n AC 
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of each of the measured signals S A .(t), S^ft) , and S AC (t) and 
can be used to estimate approximations to the desired signal" 
portions Y,.(t) , Y^t) , and Y ic (t) for all three measured 
signals S A .(t), S^t) , and S Ae (t). Using the ratiometric 
5 method, estimation of the desired signal portions Y^(t) and 
*Ac(t) of two measured signals S x .(t) and S xc (t), chosen 
correctly, is usually satisfactory to determine most 
physiological data. 

A joint process estimator 60 having two regression 

10 filters 80a and 80b is shown in Figure 9. A first regression 
filter 80a accepts a measured signal S A .(t). A second 
regression filter 80b accepts a measured signal S^t) or 
s Ac(t), depending whether the constant saturation method or 
the ratiometric method is used to determine the noise 

15 reference signal n'(t). The first and second regression 
filters 80a and 80b are independent. The backward prediction 
error b B (t) is input to each regression filter 80a and 80b, 
the input for the second regression filter 80b bypassing the 
first regression filter 80a. 

20 The second regression filter 80b comprises registers 98, 

and summing elements 108 arranged similarly to those in the 
first regression filter 80a. The second regression filter 
80b operates via an additional intermediate variable in 
conjunction with those defined by equations (48) through 

25 (58), i.e.: 

= ^.,*b(t-D + {b.(t)e* BAb (t)/ T .(t)); or (59) 
Pm,x^) - V., xe (t-1) + {b.(t)e*. fAe (t)/ 7- (t)>; and (60) 
'o.Ab(°> * °? °r (61) 
t>o,xc(°) m 0. (62) 
30 The second regression filter 80b has an error signal value 
defined similar to the first regression filter error signal 
values, e^ 1ia (t), i.e.: 

^i.AbW - e-.xbft) ~ *V*b< t > b .< t > ? or («) 
e-Mcft) = e . (AC ( t ) " «Y*b< t > b .< t >'' and (64) 
35 Zo.xbW = s >b(t) for t * 0; or (65) 
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e o. Ac ( t > = s Ac(t) for t * o. (66) 
The second regression filter has a regression coefficient' 
"..Abft) register 98 value defined similarly to the first 
regression filter error signal values, i.e.: 

- (P. fAb (t)/^(t)}; or (67) 
= (P^lt)//)^)}; (68) 
These values are used in conjunction with those intermediate 
variable values, signal values, register and register values 
defined in equations (40) through (58). These signals are 
calculated in an order defined by placing the additional 
signals immediately adjacent a similar signal for the 
wavelength Aa. 

For the ratiometric method, S AC (t) is input to the second 
regression filter 80b. The output of the second regression 
filter 80b is then a good approximation to the desired signal 
Y\c(t). For the constant saturation method, S Ab (t) is input 
to the second regression filter 80b. The output is then a 
good approximation to the desired signal Y'-^t). 

The addition of the second regression filter 80b does 
20 not substantially change the computer program subroutine 
represented by the flowchart of _ Figure 8 . Instead of an 
order update of the m th stage of only one regression filter, 
an order update of the m th stage of both regression filters 
80a and 80b is performed. This is characterized by the 
25 plural designation in the "ORDER UPDATE OF m th STAGE OF 
REGRESSION FILTER(S) " box 180 in Figure 8. Since the 
regression filters 80a and 80b operate independently, 
independent calculations can be performed in the reference 
processor and joint process estimator 60 adaptive noise 
30 canceler subroutine modeled by the flowchart of Figure 8. 

CALCULATION OF SATURATTOfrJ 
Once good approximations to the desired signals, Y' A8 (t) 
anc * Y' Ac (t) for the ratiometric method and Y' Aa (t) and Y' Ab (t) 
for the constant saturation method, have been determined by 
35 the joint process estimator 60, the saturation of ^ in a 
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volume containing A, and A^, for example, may be calculated 
according to various known methods. Mathematically, the- 
approximations to the desired signals can be written: 

Y*A.(t) « ^.05X5 6 (t) + e 6 ,,.c 6 X5 6 (t); and (69) 
Y\ e (t) « c 5 , c ^ 6 (t) + e 4<4e C 4 X, #6 (t). (70) 

for the ratiometric method using wavelengths Aa and Ac. For 
the constant saturation method, the approximations- to the 
desired signals can be written in terms of Aa and Ab as: 

" «s.jWs.«<t) + «« fl .c 4 x^ 4 <t); and (71) 
Y'^(t) « e^c^Ct) + c^c^it) . (72) 
This is equivalent to two equations having three unknowns, 
namely ^(t), c 6 (t) and 3% 6 (t) . In both the ratiometric and 
the constant saturation cases, the saturation can be 
determined by acquiring approximations to the desired signal 
portions at two different, yet proximate times t t and over 
which the saturation of A, in the volume containing A, and A 6 
does not change substantially. For example, for the desired 
signals estimated by the ratiometric method, at times t, and 



Y' A .(t,) « ^,05X5.4 (t,) + c^cyt^Ct,) (73) 

Y'AcCt,) « *5.*c c 5*S.s( t 1> + '6.^4X5.6 (74) 

Y'^t,) = ^5 ^X5 6 (t 2 ) + e^Xj^Ctj) (75) 

Y'Acd^) - «, fle <VC5 f6 (t I ) + e 4fJke ^ 6 (tj) (76) 

Then, difference signals may be determined which relate the 
signals of equation (73) through (76), i.e.: 

AY Aa = Y'^t^-Y'^tj) « ej^CjAx + c 6#x .c 6 ax; and (77) 
AYac = Y' Ac (t 1 )-Y' Ae (t 2 ) « e^C^AX + c 6 xt C 6 tX; (78) 
where ax - 6 (t,) - x^tj). The average saturation at time 
t=(t,+t 2 )/2 is: 

Saturation(t) « Cs(t)/[C5(t) + c 4 (t)) (79) 
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It will be understood that the ax term drops out from the 
saturation calculation because' of the division. Thus, 
knowledge of the thickness of the desired constituents is not 
5 required to calculate saturation. 

PULSE OXIMETRY MEASUREMENTS 

A specific example of a physiological monitor utilizing 
a processor of the present invention to determine a noise 
reference signal n'(t) for input to an adaptive noise 

10 canceler that removes erratic motion-induced undesired signal 
portions is a pulse oximeter. A pulse oximeter typically 
causes energy to propagate through a medium where blood flows 
close to the surface for example, an ear lobe, or a digit 
such as a finger, or a forehead. An attenuated signal is 

15 measured after propagation through or reflection from the 
medium. The pulse oximeter estimates the saturation of 
oxygenated blood available to the body for use. 

Freshly oxygenated blood is pumped at high pressure from 
the heart into the arteries for use by the body. The volume 

20 of blood in the arteries varies with the heartbeat, giving 
rise to a variation in absorption of energy at the rate of 
the heartbeat, or the pulse. 

Oxygen depleted, or deoxygenated, blood is returned to 
the heart by the veins along with unused oxygenated blood. 

25 The volume of blood in the veins varies with the rate of 
breathing, which is typically much slower than the heartbeat. 
Thus, when there is no motion induced variation in the 
thickness of the veins, venous blood causes a low frequency 
variation in absorption of energy. When there is motion 

30 induced variation in the thickness of the veins, the low 
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frequency variation in absorption is coupled with the erratic 
variation in absorption due to motion artifact. 

In absorption measurements using the transmission of 
energy through a medium, two light emitting diodes (LED 1 s) 
are positioned on one side of a portion of the body where 
blood flows close to the surface, such as a finger, and a 
photodetector is positioned on the opposite side of the 
finger. Typically, in pulse oximetry measurements, one LED 
emits a visible wavelength, preferably red, and the other LED 
emits an infrared wavelength. However, one skilled in the 
art will realize that other wavelength combinations could be 
used. 

The finger comprises skin, tissue, muscle, both arterial 
blood and venous blood, fat, etc., each of which absorbs 
light energy differently due to different absorption 
coefficients, different concentrations, and different 
thicknesses. When the patient is not moving, absorption is 
substantially constant except for the flow of blood. This 
constant attenuation can be determined and subtracted from 
the signal via traditional filtering techniques. When the 
patient moves, the absorption becomes erratic. Erratic 
motion induced noise typically cannot be predetermined and 
subtracted from the measured signal via traditional filtering 
techniques. Thus, determining the saturation of oxygenated 
arterial blood becomes more difficult. 

A schematic of a physiological monitor for pulse 
oximetry is shown in Figure 10. Two LED's 300 and 302, one 
LED 300 emitting red wavelengths and another LED 302 emitting 
infrared wavelengths, are placed adjacent a finger 310. A 
photodetector 320, which produces an electrical signal 
corresponding to the attenuated visible and infrared light 
energy signals is located opposite the LED^ 300 and 302. 
The photodetector 320 is connected to a single channel of 
common processing circuitry including an amplifier 330 which 
is in turn connected to a band pass filter 340. The band 
pass filter 340 passes signal into a synchronized demodulator 
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350 which has a plurality of output channels. One output 
channel is for signals corresponding to visible wavelengths- 
and another output channel is for signals corresponding to 
infrared wavelengths. 

The output channels of the synchronized demodulator for 
signals corresponding to both the visible and infrared 
wavelengths are each connected to separate paths, each path 
comprising further processing circuitry* Each path includes 
a DC offset removal element 360 and 362, such as a 
differential amplif ier, a programmable gain amplif ier 370 and 
372 and a low pass filter 380 and 382. The output of each 
low pass filter 380 and 382 is amplified in a second 
programmable gain amplifier 390 and 392 and then input to a 
multiplexer 400. - 

The multiplexer 400 is connected to an analog-to-digital 
converter 410 which is in turn connected to a microprocessor 
420. Control lines between the microprocessor 420 and the 
multiplexer 400, the microprocessor 420 and the analog-to- 
digital converter 410, and the microprocessor 420 and each 
programmable gain amplifier 370, 372, 390, and 392 are 
formed. The microprocessor 420 has additional control lines, 
one of which leads to a display 430 "and the other of which 
leads to an LED driver 440 situated in a feedback loop with 
the two LED v s 300 and 302. 

The LED's 300 and 302 each emits energy which is 
absorbed by the finger 310 and received by the photodetector 
320. The photodetector 320 produces an electrical signal 
which corresponds to the intensity of the light energy 
striking the photodetector 320 surface. The amplifier 330 
amplifies this electrical signal for ease of processing. The 
band pass filter 340 then removes unwanted high and low 
frequencies. The synchronized demodulator 350 separates the 
electrical signal into electrical signals corresponding to 
the red and infrared light energy components. A 
predetermined reference voltage, V ref , is subtracted by the DC 
offset removal element 360 and 362 from each of the separate 
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signals to remove substantially constant absorption which 
corresponds to absorption when there is no motion induced" 
undesired signal component. Then the first programmable gain 
amplifiers 370 and 372 amplify each signal for ease of 
5 manipulation. The low pass filters 380 and 382 integrate 

each signal to remove unwanted high frequency components and 
the second programmable gain amplifiers 390 and 392 amplify 
each signal for further ease of processing. 

The multiplexer 400 acts as an analog switch between 

10 the electrical signals corresponding to the red and the 
infrared light energy, allowing first a signal corresponding 
to the red light to enter the analog-to-digital convertor 410 
and then a signal corresponding to the infrared light to 
enter the analog-to-digital convertor 410. This eliminates 

15 the need for multiple analog-to-digital convertors 410. The 
analog-to-digital convertor 410 inputs the data into the 
microprocessor 420 for calculation of a noise reference 
signal via the processing technique of the present invention 
and removal of undesired signal portions via an adaptive 

20 noise canceler. The microprocessor 420 centrally controls 
the multiplexer 400, the analog-to-digital convertor 410, and 
the first and second programmable gain amplifiers 370 and 390 
for both the red and the infrared channels. Additionally, 
the microprocessor 420 controls the intensity of the LED's 

25 302 and 304 through the LED driver 440 in a servo loop to 

keep the average intensity received at the photodetector 320 
within an appropriate range. Within the microprocessor 420, 
a noise reference signal n'(t) is calculated via either the 
constant saturation method or the ratiometric method, as 

30 described above, the ratiometric method being generally 
preferred. This signal is used in an adaptive noise canceler 
of the joint process estimator type 60, described above. 

The multiplexer 400 time multiplexes, or sequentially 
switches between, the electrical signals corresponding to the 

35 red and the infrared light energy. This allows a single 

channel to be used to detect and begin processing the 
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electrical signals. For example, the red LED 300 is 
energized first and the attenuated signal is measured at the" 
photodetector 320. An electrical signal corresponding to the 
intensity of the attenuated red light energy is passed to the 
5 common processing circuitry. The infrared LED 3 02 is 
energized next and the attenuated signal is measured at the 
photodetector 320. An electrical signal corresponding to the 
intensity of the attenuated infrared light energy is passed 
to the common processing circuitry. Then, the red LED 300 is 
10 energized again and the corresponding electrical signal is 
passed to the common processing circuitry* The sequential 
energization of LED's 300 and 302 occurs continuously while 
the pulse oximeter is operating. 

The processing circuitry is divided into distinct paths 
15 after the synchronized demodulator 350 to ease time 
constraints generated by time multiplexing. In the preferred 
embodiment of the pulse oximeter shown in Figure 10, a sample 
rate, or LED energization rate, of 1000 Hz is advantageously 
employed. Thus, electrical signals reach the synchronized 
20 demodulator 350 at a rate of 1000 Hz. Time multiplexing is 
not used in place of the separate paths due to settling time 
constraints of the low pass filters 380, 382, and 384. 

In Figure 10, a third LED 304 is shown adjacent the 
finger, located near the LED's 300 and 302. The third LED 
25 304 is used to measure a third signal S Ae (t) to be used to 
determine saturation using the ratiometric method. The third 
LED 304 is time multiplexed with the red and infrared LED's 
300 and 302. Thus, a third signal is input to the common 
processing circuitry in sequence with the signals from "the 
30 red and infrared LED's 300 and 302. After passing through 
and being processed by the operational amplifier 330, the 
band pass filter 340, and the synchronized demodulator 350, 
the third electrical signal corresponding to light energy at 
wavelength Ac is input to a separate path including a DC 
35 offset removal element 364, a first programmable gain 
amplifier 374, a low pass filter 384, and a second 
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programmable gain amplifier 394. The third signal is then 
input to the multiplexer 400. 

The dashed line connection for the third LED 304 
indicates that this third LED 304 is incorporated into the 
pulse oximeter when the ratiometric method is used; it is 
unnecessary for the constant saturation method. When the 
third LED 304 is used, the multiplexer 400 acts as an analog 
switch between all three LED 300, 302, and 304 signals, if 
the third LED 304 is utilized, feedback loops between the 
microprocessor 420 and the first and second programmable gain 
amplifier 374 and 394 in the Ac wavelength path are also 
formed. 

For pulse oximetry measurements using the ratiometric 
method, the signals (logarithm converted) transmitted through 
the finger 310 at each wavelength Aa, Ab, and Ac are: 

s xa ( t) c S lr9dI ( t) = ^c^ozx * ( t) Xa c£x A (t) 
* *Hbo2.x*c£o2X v ( t) * e^^c&c v ( t) + n XM ( t) . 

(81) 

Sjj, ( t) =S lzod2 ( t) = e^ ^c^oiX* ( t) i-e^, xi> c BbX A ( t) 
+ e m>o2.xh c Zo^ v (t) * e abiXb c^x v (t) + n^it) . 

(82) 

s x c ( t) «5 iZR ( t) = e^^c^x* ( t) +e a> , lc c J £x* ( t) 
* e abo2. tcCno* v ( t) + e^, Xc cZx v * n le ( t) . 

(83) 

In equations (81) through (83), x*(t) is the lump-sum 
thickness of the arterial blood in the finger; x v (t) is the 
lump-sum thickness of venous blood in the finger; « Hb02 ia 

e Hb02.xb' c Hb0Z,xe' £ Mb,ia' c Hb,Ab' and £ Hb.Ae are absorption 

coefficients of the oxygenated and non-oxygenated hemoglobin, 
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at each wavelength measured; and c Kb02 (t) and o^{t) with the 
superscript designations A and V are the concentrations of- 
the oxygenated and non-oxygenated arterial blood and venous 
blood , respectively . 

For the ratiometric method , the wavelengths chosen are 
typically two in the visible red range, i.e., Aa and Ab, and 
one in the infrared range, i.e., Ac. As described above, the 
measurement wavelengths Aa and Ab are advantageously chosen 
to satisfy a proportionality relationship which removes the 
desired signal portion Yju(t) and Yxb(t) , yielding a noise 
reference signal n'(t). In the preferred embodiment, the 
ratiometric method is used to determine the noise reference 
signal n 1 (t) by picking two wavelengths that cause the 
desired portions Y Aa (t) and Yxbft) of the measured signals 
S*.(t) and S^t) to become linearly dependent similarly to 
equation (22); i.e. wavelengths Aa and Ab which satisfy: 

Typical wavelength values chosen are Aa « 650nm and Ab = 
685nm. Additionally a typical wavelength value for Ac is Ac 
= 940nm. By picking wavelengths Aa and Ab to satisfy 
equation (84) the venous portion of the measured signal is 
also caused to become linearly dependent even though it is 
not a portion of the desired signal. Thus, the venous 
portion of the signal is removed with the desired portion. 
The proportionality relationship between equations (81) and 
(82) which allows determination of a non-zero noise reference 
signal n 1 (t) , similarly to equation (25) is: 

n A *(t) * w^ft) • (86) 
In pulse oximetry, both equations (85) and (86) can typically 
be satisfied simultaneously. 

Figure 11 is a graph of the absorption coefficients of 
oxygenated and deoxygenated hemoglobin (c^^ and c^) vs. 
wavelength (A) . Figure 12 is a graph of the ratio of 
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absorption coefficients vs. wavelength, i.e., < H fc/ £ Hb02 vs - A 
over the range of wavelength within circle 13 in Figure 11/ 
Anywhere a horizontal line touches the curve of Figure 12 
twice, as does line 400, the condition of equation (84) is 
satisfied. Figure 13 shows an exploded view of the area of 
Figure 11 within the circle 13. Values of c Hbaz and at the 
wavelengths where a horizontal line touches the curve or 
Figure 12 twice can then be determined from the data in 
Tigure 13 to solve for the proportionality relationship of 
equation (85) . 

A special case of the ratiometric method is when the 
absorption coefficients c Hb02 and c Hb are equal at a wavelength. 
Arrow 410 in Figure 11 indicates one such location, called an 
isobestic point. Figure 13 shows an exploded view of the 
isobestic point. To use isobestic points with the 
ratiometric method, two wavelengths at isobestic points are 
determined to satisfy equation (84). 

Multiplying equation (82) by w r4 and then subtracting 
equation (82) from equation (81) , a non-zero noise reference 
signal n^t) is determined by: 

n'(t) - S Aa (t) - w^S^'t) «= n Aa (t) - w^n xb . (87) 
This noise reference signal n'(t) has spectral content 
corresponding to the erratic, motion-induced noise. When it 
is input to an adaptive noise canceler, with either the 
signals S Aa (t) and S AC (t) or S xb (t) and S AC (t) input to two 
regression filters 80a and 80b, the adaptive noise canceler 
will function much like an adaptive multiple notch filter and 
remove frequency components present in both the noise 
reference signal n»(t) and the measured signals from the 
measured signals S Aa (t) and S AC (t) or S Ab (t) and S AC (t) . Thus, 
the adaptive noise canceler is able to remove erratic noise 
caused in the venous portion of the measured signals S Aa (t), 
S A b(t), and S AC (t) even though the venous portion of the 
measured signals S Aa (t) and S Ab (t) was not incorporated in the 
noise reference signal n'(t). However, the low frequency 



absorption caused by venous blood moving through the veins is 
generally not one of the frequencies incorporated into the, 
noise reference signal n 1 (t) . Thus, the adaptive noise 
canceler generally will not remove this portion of the 
undesired signal. However, a band pass filter applied to the 
approximations to the desired signals Y» Aa (t) and Y' Ac (t) or 
Y'ib(t) and Y' ie (t) can remove this portion of the undesired 
signal corresponding to the low frequency venous absorption. 

For pulse oximetry measurements using the constant 
saturation method, the signals (logarithm converted) 
transmitted through the finger 310 at each wavelength Aa and 
Ab are: 

S Xa ( t) "S Xz9d2 ( t) «= entn^c&osx* ( t) +e Hb , Xm c£x* ( t) 
* e xba2.i*Ca>o2X v < t) ♦ *u>.i* c £x v < t) + n u ( t) . 

(88) 

s u> < c > = s iix < t) - e aM .u,c£e a x A ( t) ^c^x * ( t) 
+ ^Bbo2.xb c Zo2X v ( t) + e^, ^cZx v ( t) + n^lt). 

(89) 

For the constant saturation method, the wavelengths chosen 
are typically one in the visible red range, i.e., Aa, and one 
in the infrared range, i.e., Ab. Typical wavelength values 
chosen are Aa = 660nm and Ab «=. 940nm. Using the constant 
saturation method, it is assumed that C Kb02 (t)/C Hb (t) = 
constant. The saturation of oxygenated arterial blood 
changes slowly, if. at all, with respect to the sample rate, 
making this a valid assumption. The proportionality factor 
between equation (88) and (89) can then be written as: 

W^(t) = £ Hb02,A« C HbQ2 X ( t ) + c Hb.A» C Hb X ( t ) (go) 
c Mb02,Xb C Hb02 X ( t ) + 'Hb.AbSb*^) 

= Y' a« (t)/Y\ b ( t) ; where (91) 
n Ae (t) " « rt (t)n Ab (t) . (92) 



In pulse oximetry, it is typically the case that both 
equations (91) and (92) can be satisfied simultaneously. 

Multiplying equation (89) by t^(t) and then subtracting 
equation (89) from equation (88), a non-zero noise reference 
signal n'(t) is determined by: 

n'(t) ■= S At (t) - u^(t)S Ab (t) (93) 

* e HbQ2. x*Cmx>*x v ( t) ♦e^ Xm c£x v ( t) +n XM ( t) 

(94) 

The constant saturation assumption does not cause the 
venous contribution to the absorption to be canceled along 
with the desired signal portions Y Aa (t) and Y xb (t) , as did the 
relationship of equation (84) used in the ratiometric method. 
Thus, frequencies associated with both the low frequency 
modulated absorption due to venous absorption when the 
patient is still and the erratically modulated absorption due 
to venous absorption when the patient is moving are 
represented in the noise reference signal n'(t). Thus, the 
adaptive canceler can remove both erratically modulated 
absorption due "to venous blood in the finger under motion and 
the constant low frequency cyclic absorption of venous blood* 

Using either method, a noise reference signal n*(t) is 
determined by the processor of the present invention for use 
in an adaptive noise canceler which is defined by software in 
the microprocessor. The preferred adaptive noise canceler is 
the joint process estimator 60 described above. 

Illustrating the operation of the ratiometric method of 
the present invention , Figures 14 , 15 and 16 show signals 
measured for use in determining the saturation of oxygenated 
arterial blood using a reference processor of the present 
invention which employs the ratiometric method, i.e., the 
signals S Aa (t) « S Ared1 (t), S xb (t) « S Arcd2 (t) , and S Ac (t) = 
S A u(t). A first, segment . 14a, 15a, and 16a of each of the 
signals is relatively undisturbed by motion artifact, i.e., 

-55- 



# • 

the patient did not move substantially during the time period 
in which these segments were measured. These segments 14a,* 
15a, and 16a are thus generally representative of the desired 
plethysmography waveform at each of the measured 
wavelengths. A second segment 14b, 15b, and 16b of each of 
the signals is affected by motion artifact, i.e. , the patient 
did move during the time period in which these segments were 
measured. Each of these segments 14b, 15b, and 16b shows 
'large motion induced excursions in the measured signal. A 
third segment 14c, 15c, and 16c of each of the signals is 
again relatively unaffected by motion artifact and is thus 
generally representative of the desired plethysmography 
waveform at each of the measured wavelengths. 

Figure 17 shows the noise reference signal n 1 (t) « n Aa 
- w r 4*Ub(t) , as determined by a reference processor of the 
present invention utilizing the ratiometric method. As 
discussed previously, the noise reference signal n»(t) is 
correlated to the undesired signal portions n Aa , n^, and n Ac . 
Thus, a first segment 17a of the noise reference signal n» (t) 
is generally flat, corresponding to the fact that there is 
very little motion induced noise in the first segments 14a, 
15a, and 16a of each signal. A second segment 17b of the 
noise reference signal n'(t) exhibits large excursions, 
corresponding to the large motion induced excursions in each 
of the measured signals. A third segment 17c of the noise 
reference signal n 1 (t) is generally flat, again corresponding 
to the lack of motion artifact in the third segments 14a, 
14b, and 14c of each measured signal. 

Figures 18 and 19 show the approximations Y' Aa (t) and 
Y* Ac (t) to the desired signals Y Aa (t) and Y AC (t) as estimated 
by the joint process estimator 60 using a noise reference 
signal n* (t) determined by the ratiometric method. Note that 
the scale of Figures 14 through 19 is not the same for each 
figure to better illustrate changes in each signal. Figures 
18 and 19 illustrate the effect of the joint process 
estimator adaptive noise canceler using the noise reference 



signal n*(t) as determined by the reference processor of the 
present invention using the ratiometric method. Segments 18b 
and 19b are not dominated by motion induced noise as were 
segments 14b, 15b, and 16b of the measured signals. 
Additionally, segments 18a, 19a, 18c, and 19c have not been 
substantially changed from the measured signal segments 14a, 
15a, 16a, 14c, 15c, and 16c vhere there was no motion induced 
noise. 

Illustrating the operation of the constant saturation 
method of the present invention, Figures 20 and 21 show 
signals measured for input to a reference processor of the 
present invention which employs the constant saturation 
method, i.e., the signals S Aa (t) = S Arod (t) and S^ft) = S AIR (t). 
A first segment 20a and 21a of each of the signals is 
relatively undisturbed by motion artifact, i.e., the patient 
did not move substantially during the time period in which 
these segments were measured. These segments 20a and 21a are 
thus generally representative of the desired plethysmography 
waveform at each of the measured wavelengths. A second 
segment 20b and 21b of each of the signals is affected by 
motion artifact, i.e., the patient did move during the time 
period in which these segments were measured. Each of these 
segments 20b and 21b shows large motion induced excursions in 
the measured signal. A third segment 20c and 21c of each of 
the signals is again relatively unaffected by motion artifact 
and is thus generally representative of the desired 
plethysmographic waveform at each of the measured 
wavelengths . 

Figure 22 shows the noise reference signal n f (t) *= n A ,(t) 
- w i4nAb(t) , as determined by a reference processor of the 
present invention utilizing the constant saturation method. 
Again, the noise reference signal n 1 (t) is correlated to the 
undesired signal portions n Aa and n^. Thus, a first segment 
22a of the noise reference signal n f (t) is generally flat, 
corresponding to the fact that there is very little motion 
induced noise in the first segments 20a and 21a of each 
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signal. A second segment 22b of the noise reference signal 
n'(t) exhibits large excursions, corresponding to the large" 
motion induced excursions in each of the measured signals. 
A third segment 22c of the noise reference signal n'(t) is 
generally flat, again corresponding to the lack of motion 
artifact in the third segments 20b and 21c of each measured 
signal. 

Figures 23 and 24 show the approximations Y' Aa (t) and 
^■^(t) to the desired signals Y w (t) and Y Ab (t) as estimated 
by the joint process estimator 60 using a noise reference 
signal n'(t) determined by the constant saturation method. 
Note that the scale of Figures 20 through 24 is not the same 
for each figure to better illustrate changes in each signal. 
Figures 23 and 24 illustrate the effect of the joint process 
estimator adaptive noise canceler using the noise reference 
signal n^t) as determined by a reference processor of the 
present invention utilizing the constant saturation method. 
Segments 23b and 24b are not dominated by motion induced 
noise as were segments 20b and 21b of the measured signals. 
Additionally, segments 23a, 24a, 23c, and 24c have not been 
substantially changed from the measured signal segments 20a, 
21a, 20c, and 21c where there was no motion induced noise. 
METHOD FOR ESTIMATING DESIRED PORTIONS OF 
MEASURED SIGNALS IN A PULSE OXIMETER 

A copy of a computer program subroutine, written in the 
C programming language, calculates a noise reference signal 
n 1 (t) using the ratiometric method and, using a joint process 
estimator 60, estimates the desired signal portions of two 
measured signals, each having an undesired portion which is 
correlated to the noise reference signal n'(t) and one of 
which was not used to calculate the noise reference signal 
n*(t), is appended in Appendix A. For example, S A- (t) « 
S A redi(t) «= S^ 50rw (t) and S AC (t) «= S AIR (t) - S xWkm (t) can be input 
to the computer subroutine. One skilled in the art will 
realize that S Att (t) - S Arcd2 (t) - S A685fW (t) and S AC (t) = S AlR (t) = 
S A 94o™(t) will also work. This subroutine is one way to 
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implement the steps illustrated in the flowchart of Figure 8 
for a monitor particularly adapted for pulse oximetry. 

The program estimates the desired signal portions of two 
light energy signals, one preferably corresponding to light 
5 in the visible red range and the other preferably 

corresponding to light in the infrared range such that a 
determination of the amount of oxygen available to the body, 
or the saturation of oxygen in the arterial blood, may be 
•made. The calculation of the saturation is performed in a 
10 separate subroutine. Various methods for calculation of the 
oxygen saturation are known to those skilled in the art. One 
such calculation is described in the articles by G.A. Hook, 
et al, and Michael R. Neuman cited above. Once the 
concentration of oxygenated hemoglobin and deoxygenated 
15 hemoglobin are determined, the value of the saturation is 
determined similarly to equations (73) through (80) wherein 
measurements at times t, and are made at different, yet 
proximate times over which the saturation is relatively 
constant. For pulse oximetry, the average saturation at time 
t = (^+1^/2 is then determined by: 

Saturation(t) = C^it) / {C Hb02 (t) + C Hb (t)]. (95) 

e Ht. Lm~ e iO>. JU> ( a ) 
- al Xb 



20 



A y Aa 

e ab.Xa~ e Hb02.Xa~^ e a>.Xb~ e Sb02.Xb^ ( A y * ) 



(96) 

Using the ratiometric method, three signals S A ,(t) , 
Svb(t) , and S Ae (t) are input to the subroutine. S Aa (t) and 
25 S ib (t) are used to calculate the noise reference signal n' (t) . 
As described above, the wavelengths of light at which S A ,(t) 
and S Ab (t) are measured are chosen to satisfy the relationship 
of equation (84). Once the noise reference signal n«(t) is 
determined, the desired signal portions Y Aa (t) and Y ic (t) of 
the measured signals S Aa (t) and S ic (t) are estimated for use 
in calculation of the oxygen saturation. 



30 
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The correspondence of the program variables to the 
variables defined in the discussion of the joint process * 
estimator is as follows: 

A^(t) = nc[] .Delta 
5 r f,^t> = nc[].fref 

r tmW = nc[].bre£ 
f«(t) = nc[],ferr 
b a (t) = nc[].berr 
^(t) = nc[].Fswsqr 
10 £ B (t) = nc[].Bswsqr 

7 a (t) ■= nc[]. Gamma 

" nc[].Roh_a 
'■xcf*) = nc[].Roh_c 
e « # A« (t) " nc tl- err _ a 
15 e «\c(*) " nc[].err_c 

V\a( t ) " nc[].K_a 
Vxcft) - nc[].K_c 

A first portion of the program performs the 

20 initialization of the registers 90 , 92 , 96, and 98 and 

intermediate variable values as in the "INITIALIZE NOISE 
CANCELER" box 120 and equations (40) through (44) and 
equations (61) , (62) , (65) / and (66). A second portion of 
the program performs the time updates of the delay element 

25 variables 110 where the value at the input of each delay 
element - variable 110 is stored in the delay element variable 
110 as in the "TIME UPDATE OF [Z" 1 ] ELEMENTS" box 130, 

A third portion of the program calculates the noise 
reference signal, as in the "CALCULATE NOISE REFERENCE 

30 (n'(t)) FOR TWO MEASURED SIGNAL SAMPLES" box 140 using the 

proportionality constant w r4 determined by the ratiometric 
method as in equation (85) . 

A fourth portion of the program performs the zero-stage 
update as in the H ZERO-STAGE UPDATE* box 150 where the zero- 

35 stage forward prediction error f 0 (t) and the zero-stage 
backward prediction error b 0 (t) are set equal to the value of 
the noise reference signal n* (t) just calculated. 
Additionally, zero-stage values of intermediate variables 
^ 0 (t) and 0 o (t) (nc[]*Fswsqr and nc[].Bswsqr in the program) 

40 are calculated for use in setting register 90, 92, 96, and 98 
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values in the least-squares lattice predictor 70 and the 
regression filters 80a and 80b. 

A fifth portion of the program is an iterative loop 
wherein the loop counter, m, is reset to zero with a maximum 
5 of m « NC_CELLS, as in the "m=0" box 160 in Figure 8. 

NC_CELLS is a predetermined maximum value of iterations for 
the loop. A typical value of NC_CELLS is between 60 and 80, 
for example* The conditions of the loop are set such that 
the loop iterates a minimum of five times and continues to 

10 iterate until a test for conversion is met or m « NC_CELLS. 

The test for conversion is whether or not the sum of the 
weighted sum of forward prediction errors plus the weighted 
sum of backward prediction errors is less than a small 
number, typically 0.00001 (i.e, f m (t)+0 m (t) s 0.00001). 

15 A sixth portion of the program calculates the forward 

and backward reflection coefficient r m f (t) and b (t) register 
90 and 92 values (nc[].fref and nc[].bref in the program) as 
in the "ORDER UPDATE m th -STAGE OF LSL- PREDICTOR" box 170 and 
equations (49) and (50) . Then forward and backward 

20 prediction errors f m (t) and b B (t) (nc[].ferr and nc[].berr in 
the program) are calculated as in equations (51) and (52) . 
Additionally, intermediate variables ^(t) , fi m (t) and 7 w (t) 
(nc[] .Fswsqr, nc[].Bswsqr, nc[]. Gamma in the program) are 
calculated, as inequations (53), (54), and (55). The first 

25 cycle of the loop uses the values for nc[0]. Fswsqr and 
nc[0].Bswsqr calculated in the ZERO-STAGE UPDATE portion of 
the program. 

A seventh portion of the program, still within the loop, 
calculates the regression coefficient * B#Aa (t) and * -#xc (t) 

30 register 96 and 98 values (nc[].K_a and nc[].K_c in the 
program) in both regression filters, as in the "ORDER UPDATE 
m th STAGE OF REGRESSION FILTER (S) " box 180 and equations (57) 
through (68). Intermediate error signals and variables 
'Wt* ' ' Pm.*™ > and Vxc^) <nc[ J .erra and 

35 nc[].err_c, nc[].roh_a, and nc[ ] .roh_c in the subroutine) are 



-61- 




also calculated as in equations (58), (64), (56), and (60), 
respectively* 

The test for convergence of the joint process estimator 
is performed each time the loop iterates, analogously to the 
5 "DONE" box 190. If the sum of the weighted sums of the 

forward and backward prediction errors f m (t)+p m (t) is less 
than or equal to 0.00001, the loop terminates. Otherwise, 
,the sixth and seventh portions of the program repeat. 

When either the convergence test is passed or m «= 

10 NC_CELLS, an eighth portion of the program calculates the 
output of the joint process estimator 60 adaptive noise 
canceler as in the "CALCULATE OUTPUT" box 200. This output 
is a good approximation to both of the desired signals Y , At (t) 
and Y , AC (t) for the set of samples S Aa (t), S Ab (t) , and S Ac (t) 

15 input to the program. After many sets of samples are 

processed by the joint process estimator, a compilation of 
the outputs provides output waves which are good 
approximations to the plethysmography wave at each 
wavelength, Aa and Ac. 

20 Another copy of a computer program subroutine, written 

in the C programming language, which calculates a noise 
reference signal n 1 (t) using the constant saturation method 
and, using a joint process estimator 60, estimates a good 
approximation to the desired signal portions of two measured 

25 signals, each having an unde sired portion which is correlated 
to the noise reference signal n'(t) and each having been used 
to calculate the noise reference signal n' (t) , is appended in 
Appendix B. This subroutine is another way to implement the 
steps illustrated in the flowchart of Figure 8 for a monitor 

30 particularly adapted for pulse oximetry. The two signals are 
measured at two different wavelengths Aa and Ab, where Aa is 
typically in the visible region and Ab is typically in the 
infrared region. For example, in one embodiment of the 
present invention, tailored specifically to perform pulse 
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oximetry using the constant saturation method , Aa = 660nm and 
Xh = 94 0nnu 

The correspondence of the program variables to the 

variables defined in the discussion of the joint process 

5 estimator is as follows: 

A.(t) = nc[] .Delta 
r f .(t) - nc[].fref 

^'.(t) " nc[].bref 

f.'(t) - nc[].ferr 
10 b.(t) - nc[].berr 

*^(t) = nc[].Fsvsgr 

p m (t) « nc[).Bsvsgr 

7(t) « nc(] .Gamma 

*m x«( t ) " nc[].Roh_a 
15 P ^ b (t) «= nc[].Rohb 

e . # Aa( t ) " nc[].err_a 

e . xb( t ) * nc[].err_b 
(t > " nc[].X_a 
- nc[].Kb 

20 

First and second portions of the subroutine are the same 
as the first and second portions of the above described 
subroutine tailored for the ratiometric method of determining 
the noise reference signal n 1 (t) . 

25 A third portion of the subroutine calculates the noise 

reference signal , as in the "CALCULATE NOISE REFERENCE 
(n»(t)) FOR TWO MEASURED SIGNAL SAMPLES" box 140 for the 
signals S Aa (t) and S xb (t) using the a proportionality constant 
^(t) determined by the constant saturation method as in 

30 equations (90) and (91). The saturation is calculated in a 
separate subroutine and a value of w^(t) is Imported to the 
present subroutine for estimating the desired portions Y Xa (t) 
and Y Ab (t) of the composite measured signals S Aa (t) and S Ab (t) . 
Fourth, fifth, and sixth portions of the subroutine are 

35 similar to the fourth, fifth, and sixth portions of the above 
described program tailored for the ratiometric method. 
However, the signals being used to estimate the desired 
signal portions Y Aa (t) and Y^t) in the present subroutine 
tailored for the constant saturation method, are S Aa (t) and 



-63- 



SibttJr the same signals that were used to calculate the noise 
reference signal n^t). 

A seventh portion of the program, still within the loop 
begun in the fifth portion of the program, calculates the 
regression coefficient register 96 and 98 values <c BX| (t) and 
* -#A b(t) (nc[].K_a and nc[].K_b in the program) in both 
regression filters* as in the "ORDER UPDATE m th STAGE OF 
REGRESSION FILTER(S)" box 180 and equations (57) through 
(67). Intermediate error signals and variables Aa (t), 
^xbf*)' and ^Ab^) (nc[].err_a and nc[].err_b, 

nc [ ] . roh_a , and nc [ ] . roh_b in the subroutine) are also 
calculated as in equations (58) , (63), (56) , and (59) , 
respectively. 

The loop iterates until the test for convergence is 
passed, the test being the same as described above for the 
subroutine tailored for the ratiometric method. The output 
of the present subroutine is a good approximation to the 
desired signals Y' Aa (t) and Y'^ft) for the set of samples 
S A «(t) and S xb (t) input to the program. After approximations 
to the desired signal portions of many sets of measured 
signal samples are estimated by the joint process estimator, 
a compilation of the outputs provides waves which are good 
approximations to the plethysmography wave at each 
wavelength, Aa and Xb. The estimating process of the 
iterative loop is the same in either subroutine, only the 
sample values S Aa (t) and S AC (t) or S Aa (t) and S Ab (t) input to 
the subroutine for use in estimation of the desired signal 
pprtions Y Aa (t) and Y Ac (t) or Y Aa (t) and Y^Ct) and how the 
noise reference signal n»(t) is calculated are different for 
the ratiometric method and the constant saturation methods. 

Independent of the method used, ratiometric or constant 
saturation, the approximations to the desired signal values 
Y' Aa (t) and Y' Ac (t) or Y' Aa (t) and Y f Ab (t) are input to a 
separate subroutine in which the saturation of oxygen in the 
arterial blood is calculated. If the constant saturation 



method is used, the saturation calculation subroutine also 
determines a value for the proportionality constant t^(t) as* 
defined in equations (90) and (91) and discussed above. The 
concentration of oxygenated arterial blood can be found from 
the approximations to the desired signal values since the 
desired signals are made up of terms comprising x(t), the 
thickness of arterial blood in- the finger; absorption 
coefficients of oxygenated and de-oxygenated hemoglobin , at 
each measured wavelength; and C^^Jt) and C^ft) , the 
concentrations of oxygenated and de-oxygenated hemoglobin, 
respectively. The saturation is a ratio of the concentration 
of one constituent, Aj, with respect to the total 
concentration of constituents in the volume containing and 
A 6 . Thus, the thickness, x(t), is divided out of the 
saturation calculation and need not be predetermined. 
Additionally, the absorption coefficients are constant at 
each wavelength. The saturation of oxygenated arterial blood 
is then determined as in equations (95) and (96) . 

While one embodiment of a physiological monitor 
incorporating a processor of the present invention for 
determining a noise reference signal for use in an adaptive 
noise canceler to remove erratic noise components from a 
physiological measurement has been described in the form of 
a pulse oximeter, it will be obvious to one skilled in the 
art that other types of physiological monitors may also 
employ the above described techniques for noise reduction on 
a composite measured signal in the presence of noise. 

Furthermore, it will be understood that transformations 
of measured signals other than logarithmic conversion and 
determination of a proportionality factor which allows 
removal of the desired signal portions for determination of 
a noise reference signal are possible. Additionally, 
although the proportionality factor w has been described 
herein as a ratio of a portion of a first signal to a portion 
of a second signal, a similar proportionality constant 
determined as a ratio of a portion of a second signal to a 
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portion of a first signal could equally well be utilized in 
the processor of the present invention. In the latter case," 
a noise reference signal would generally resemble n' (t) = 
n Ab (t) - «n Aa (t) . 

It will also be obvious to one skilled in the art that 
for most physiological measurements, two wavelengths may be 
determined which will enable a signal to be measured which is 
indicative of a quantity of a component about which 
information is desired. Information about a constituent of 
any energy absorbing physiological material may be determined 
by a physiological monitor incorporating a signal processor 
of the present invention and an adaptive noise canceler by 
determining wavelengths which are absorbed primarily by the 
constituent of interest. For most physiological 

measurements, this is a simple determination. 

Moreover, one skilled in the art will realize that any 
portion of a patient or a material derived from a patient may 
be used to take measurements for a physiological monitor 
incorporating a processor of the present invention and an 
adaptive noise canceler. Such areas include a digit such as 
a finger, but are not limited to a finger. 

One skilled in the art will realize that many different 
types of physiological monitors may employ a signal processor 
of the present invention in conjunction with an adaptive 
noise canceler. Other types of physiological monitors 
include, but are in not limited to, electron cardiographs, 
blood pressure monitors, blood gas saturation (other than 
oxygen saturation) monitors, capnographs, heart rate 
monitors, respiration monitors, or depth of anesthesia 
monitors. Additionally, monitors which measure the pressure 
and quantity of a substance within the body such as a 
breathalizer, a drug monitor, a cholesterol monitor, a 
glucose monitor, a carbon dioxide monitor, a glucose monitor, 
or a carbon monoxide monitor may also employ the above 
described techniques for removal of undesired signal 
portions. 
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Furthermore, one skilled in the art will realize that 
the above described techniques of noise removal from a* 
composite signal including noise components can also be 
performed on signals made up of reflected energy, rather than 
5 transmitted energy. One skilled in the art will also realize 

that a desired portion of a measured signal of any type of 
energy , including but not limited to sound energy, X-ray 
energy , gamma ray energy, or light energy can be estimated by 
the noise removal techniques described above. Thus, one 

10 skilled in the art will realize that the processor of the 
present invention and an adaptive noise canceler can be 
applied in such monitors as those using ultrasound where a 
signal is transmitted through a portion of the body and 
reflected back from within the body back through this portion 

15 of the body. Additionally, monitors such as echo 
cardiographs may also utilize the techniques of the present 
invention since they too rely on transmission and reflection. 

While the present invention has been described in terms 
of a physiological monitor, one skilled in the art will 

20 realize that the signal processing techniques of the present 
invention can be applied in many areas, including but not 
limited to the processing of a physiological signal. The 
present invention may be applied in any situation where a 
signal processor comprising a detector receives a first 

25 signal which includes a first desired signal portion and a 
first undesired signal portion and a second signal which 
includes a second desired signal portion and a second 
undesired signal portion. The first and second signals 
propagate through a common medium and the first and second 

30 desired signal portions are correlated with one another. 

Additionally, at least a portion of the first and second 
undesired signal portions are correlated with one another due 
to a perturbation of the medium while the first and second 
signals are propagating through the medium. The processor 

35 receives the first and second signals and combines the first 

and second signals to generate a noise reference signal in 
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which the primary component is derived from the first and 
second undesired signal portions. Thus, the signal processor^" 
of the present invention is readily applicable to numerous 
signal processing areas. 
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/************************************** *********************** ******** A#A ^^^ 
************************** APPeut*^ A ****************************** 
************************* i^east Square Lattice ************************* # ^^ # 
************************ Noise Cancelling **************************** 

/* Example for ratiometric approach to noise cancelling */ 
/define LAMBDA 0.95 

void OxiLSL_NC( int reset, 

int passes, 

int *signal_l, 

int *signal_2, 

int *signal~3 , 

int *target_l, 

int *target_2) { 

int i, ii, k, in, n, contraction; 
static int *s_a, *s_b, *s_c, *out_a, *out_c; 
static float Delta_sqr, scale, noise_ ref; 

if ( reset ■*= TRUE) { 
s_a « signal_l; 
s_b *= signal_2; 
6_c »= signal_3; 
out_a « target_l; 
out_c ■» target_2; 
factor ■= 1,5; 
scale ■= 1.0 /4160.0; 

/* noise canceller initialization at time t=0 */ 

nc[0] .berr «= 0.0; 
nc[0] •Gamma *= 1.0; 

for (o=0; m<NC_CELLS; m++) { 

nc[m].err_a «= 0.0; 

nc[m].err_b = 0.0; 

nc[m].Roh_a ■= 0.0; 

nc[m] .Roh_c ■= 0.0; 

nc[m] .Delta ■= 0.0; 

nc[m] .Fsvsgr « 0.00001; 

nc[m] .Bswsqr 0.00001; 



/* ■ END INITIALIZATION 

for(k=0; k<passes; k++) { 
contraction ■= FALSE; 

for(m=0; a< NC_CELLS; m++) { /* Update delay elements 

nc[m].berrl « nc[o].berr; 
nc[m] .Bswsqrl nc[m] .Bsvsqr; 

} 



noise_ref «= factor * log(l.O - (*s_a) * scale) 

- log(1.0 - (*s_b) * scale) ; 

nc[0].err_a «= log(1.0 - (*s_b) * scale); 

nc[0].err__b = log(1.0 - (*s_c) * scale); 
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++s_b; 
++s_c; 

nc[0] .f err 
nc(0) .berr 
nc[0] .Fswsqr 
nc[0] .Bswsqr 



noiee_ref ; 
noise~ref ; 

LAMBDA * nc[0]. Fswsqr + noise_ref * noisejref; 
nc[0) .Fswsqr; 



/* Order Update */ 

for(n«l;( n < NC_CELLS) 44 (contraction — FALSE); n++) { 
/* Adaptive Lattice Section */ * 



a - n-1; 
ii- n-l; 

nc[a] .Delta *- 
nc[m] .Delta +« 
Delta_sqr ■ 



LAMBDA; 

nc[m].berrl * nc[m].ferr / nc[a]«.Gaaaa ; 
nc [a]. Delta * nc [a]. Delta; 



nc[n].fref « -nc[aj .Delta / nc[a) .Bswsqrl; 
nc[n].bref - -nc[a]. Delta / nc[a) .Fswsqr; 

nc[n].ferr ~ nc[m].ferr + nc[n].fref * nc[m] .berrl; 
nc[n].berr « nc[a]. berrl + nc[n].bref * nc[a].ferr; 

nc[n] .Fswsgr «= nc[a]. Fswsqr - Delta_sqr / nc [a] .Bswsqrl; 
nc[n]. Bswsqr ~ nc[a) .Bswsqrl - Delta_sqr / nc[a] .Fswsqr; 

if( (nc[n]. Fswsqr + nc[n] .Bswsqr) > 0.00001 | | (n < 5) ) { 

nc[n]. Gamma « nc[m].Gaama - nc[a]. berrl * nc[a]. berrl / nc [a] .Bswsqrl; 
if (nc[n] .Gamma < 0.05) nc[n) .Gamma = 0.05; 
if (nc[nj .Gamma > 1.00) nc[n]. Gamma = 1.00; 

/* Joint Process Estimation Section */ 

nc[m].Roh_a *= LAMBDA; 

nc[m].Roh_a += nc[a].berr * nc[m].err_a / nc[a]. Gamma ; 
nc[a].k_a «= nc[a].Roh_a / nc[a] .Bswsqr; 
nc[n].err_a « nc[a].err_a - nc[a].k_a * nc[a].berr; 

nc[a].Roh_c LAMBDA; 

nc[a].Roh_c +- nc[a].berr * nc[a].err_b / nctaJ.Gaama ; 
nc[a].k_c ■= nc[a].Roh_c / nc[a] .Bswsqr; 
. nc[n].err_b « nc(a].err_b - nc[m].k_c * nc[a].berr; 

} 

else { 

contraction «= TRUE; 
for(i«n; i<NC_CELLS; i++) { 



nc[ij .err_a 
nc[i] .Roh_a 
nc[i) .err_b 
nc[i] .Roh_c 
nc[i] .Delta 
nc[i] .Fswsqr 
nc[i] .Bswsqr 
nc[i] .Bswsqrl 



0.0; 
0.0; 
0.0; 
0.0; 
0.0; 
0.00001; 
0.00001; 
0.00001; 
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> 

> 

) 

*out_a++ - (int) ( (-exp(nc[ii] .err a) +1.0) / scale) ; 

*out_c++ - (int) ( (-exp(nc[ii] .err~b) +1.0) / scale) ; 



) 

} 

/********************** Least Square Lattice *************************** 
************************* ****************************** 
*****************************************************************************/ 



41 
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/******* ******************************************************* *************** 
************************** ^ppejoi^ijc *H **************************** ## 

************************* Least Square Lattice ***************************** 
************************ Hoise Cancelling **************************** / 

/* Example for constant saturation approach to noise cancelling */ 
/define LAMBDA 0.95 

void 0xiLSL_NC( int reset, 

lnt passes, 

int sat_f actor , 

int *signal_l, 

int *signal_2, 

int * target 1, 

int *target~2) { 

int i, ii, X, a, n, contraction; 
static int *6_a, *s_ b, *out_a, *out_b; 
static float Delta_sgr, scaleT noisejref ; 

if ( reset — TRUE){ 
s_a « signal_l; 
sjb ■= signal_2; 
out_a «= target_l; 
out_b « target_2; 
scale « 1,0 /4160.0; 

/* noise canceller initialization at time t«0 */ 

nc(0].berr « 0.0; 
nc[0] .Gamma = 1.0; 

for(m=0; m<NC_CEL»L»S ; m++) { 

nc[m].err_a «= 0.0; 

nc[m].err_b -=0.0; 

^ nc[m] .Roh_a *= 0.0; 

nc[m].Roh_b =0.0; 

nc[m]. Delta « 0.0; 

nc[m] .Fswscpr ■= 0.00001; 

nc[m].Bsvsgr ■= 0.00001; 



/* END INITIALIZATION 

for(k«0; k<passes; k++) { 
contraction FALSE; 

for(m=0; m< NC_CELLS; m++) { /* Update delay elements 

nc[m].berrl ■= nc[m].berr; 
nc[m] .Bsvsgrl « nc[m] .Bsw6gr; 

} 



noise_ref « sat_factor * log(1.0 - (*s_a) * scale) 

- log (1.0 - (*s_b) * scale) ; 
nc[0].err_a = log (1.0 - (*s_a) * scale); 
nc(0].err_b log(1.0 - (*s_b) * scale) ; 

++s_a ; 
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* 




++s_b; 



nc[0].ferr «= noise_ref ; 
nc[0].berr ■= noise_ref ; 

nc[0). Fsvsqr «= LAMBDA * nc[0]. Fsvsqr + noise_ref * noise_ref; 
nc(O). Bsvsqr «= nc[0] .Fsvsqr; 



for(n-l;( n < NC_CELLS) (contraction — = FALSE); n++) { " 

/* Adaptive Lattice Section */ 

m - n-1; 
ii- n-1; 

nc[m]. Delta *= LAMBDA; 

nc[m] .Delta nc[m]. berrl * nc[m].ferr / nc[m]. Gamma ; 

Delta_sqr « nc[m] .Delta * nc[m]. Delta; 

nc[nj.fref « -nc[m] .Delta / nc[m] .Bswsqrl; 
nc[n].bref « -nc[m). Delta / nc[&] .Fsvsqr; 

nc[n].ferr «= nc[m].ferr + nc[n].fref * nc[m]. berrl; 
nc[n].berr - nc[m]. berrl + nc[n].bref * nc[m).ferr; 

nc[n]. Fsvsqr = nc[m] .Fsvsqr - Delta_sgr / nc[m] .Bsvsqrl; 
nc[n]. Bsvsqr «= nc[m) .Bsvsqrl - Delta_sqr / nc[m) .Fsvsqr; 

if( (nc[n] .Fsvsqr + nc[n] .Bsvsqr) > 0.00001 | | (n < 5) ) { 

nc[n] .Gamma ■= nc[m] .Gamma - nc[m]. berrl * nc[mj.berrl / nc[m] .Bsvsqrl; 
if (nc(n) .Gamma < 0.05) nc[n] .Gamma ■= 0.05; 
if (nc[n] .Gamma > 1.00) nc[n]. Gamma = 1.00; 

/* Joint Process Estimation Section */ 

nc[m].Roh_a *= LAMBDA; 

nc[m].Roh_a += nc[m].berr * nc[m).err_a / nc[m]. Gamma ; 
nc[m].k_a « nc[m].Roh_a / nc[m] .Bsvsqr; 
nc[n].err_a « nc[m).err_a - nc[m].k_a * nc[m].berr; 

nc[m].Roh_b *«= LAMBDA; 

nc[m].Roh_b nc[m].berr * nc[m].err_b / nc[m]. Gamma ; 
nc[m).k_b «= nc[m].Roh_b / nc[m] .Bsvsqr; 
nc[n).err_b •= nc[m).err_b - nc[m].k_b * nc[m].berr; 



else { 

contraction « TRUE; 
for(i«n; i<NC_CELLS; i++) { 

nc[i].err_a « 0.0; 

nc[i].Roh_a « 0.0; 

nc[i] .err_b « 0.0; 

nc[i].Roh_b = 0.0; 

nc£i]. Delta ■= 0.0; 

nc [ i ] . Fsvsqr ■= 0 . 00001 ; 

nc[i]. Bsvsqr «= 0.00001; 

nc[i] .Bsvsqrl = 0.00001; 



/* Order Update 



} 



} 
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} 

*out_a++ « (int) ( (-exp(nc[iij .err_a) +1.0) / scale) ; 
*out_b++ «s (int) ( (-exp(nc[ii].err_b) +1.0) / scale) ; 

} 

} 

/********************** Least Square Lattice A************************** 
************************* ***************************** 
*************«*«******«***«««*«**«***«*******«**« A *********«*******«***«« Aft#v 
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